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I .  INTRODUCTION 


In  this  paper  a  description  of  a  reverse  annular  regenerative 

liquid  propellant  gun  code  is  given.  The  basic  governing  equations  for 

the  lumped  parameter  model  are  similar  to  those  developed  for  the  in- 
1-T 

line  model,  and  the  derivations  are  not  repeated  in  detail.  Changes 
in  the  equations  are  described,  and  the  reasons  for  the  changes  are 
given.  A  complete  description  of  the  present  group  of  options  is  given. 
The  code  is  compared  to  experimental  data  for  a  30mm  gun  fixture 
developed  by  General  Electric.  &  >-•  >«.  j  C  *  7  .-*~ 

II.  THE  REGENERATIVE  LIQUID  PROPELLANT  GUN 

Figure  1  shows  a  diagram  of  a  generic  regenerative  liquid 
propellant  gun  with  a  reverse  piston.  The  monopropellant  in  the  liquid 
reservoir  is  around  the  gun  tube.  A  primer  is  ignited  and  injects  hot 
gas  into  the  combustion  chamber.  As  the  chamber  is  pressurized,  the 
piston  is  pushed  to  the  right.  Because  of  the  piston  area  differential 
between  the  two  regions,  the  piston  will  move  forward  even  when  the 
liquid  pressure  is  higher  than  the  intermediate  chamber  pressure.  The 
pressure  differential  forces  liquid  propellant  through  the  holes  in  the 
piston  face  and  through  the  holes  directly  into  the  gun  tube.  The 
propellant  ignites  and  burns  in  the  chamber  and  in  the  tube.  Hot  gas  may 
flow  from  the  Intermediate  chamber  into  the  tube  or  from  the  tube  into 
the  intermediate  chamber,  depending  on  the  relative  size  of  the 
injection  areas.  The  increasing  chamber  pressure  accelerates  the 
piston,  leading  to  more  rapid  injection.  Eventually,  the  pressure 
pushes  the  projectile  down  the  gun  tube. 

Figure  2  shows  the  design  of  the  GE  fixture.  The  booster  at  the 
left  holds  the  igniter.  Hot  gas  is  vented  from  a  series  of  holes  around 
the  booster.  There  are  large  slots  in  the  stub  case  holder  and  the  end 
of  the  short  barrel  ,  so  the  hot  gas  easily  penetrates  to  the  piston 
face.  The  piston  face  has  no  injection  vents,  which  simplifies 
the  propellant  loading  procedure.  The  liquid  propellant  is  initially 
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Figure  1.  A  Generic  Reverse  Annular  Regenerative  Liquid  Propellant  Gun, 
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Figure  2.  GE  30nun  Reverse  Annular  Regenerative  Liquid  Propellant  Gun. 


around  the  short  barrel  between  the  piston  and  the  valve  head.  The 
injection  occurs  through  a  series  of  "Y"  shaped  holes  into  the  gun  tube. 
The  projectile  initially  blocks  these  holes.  As  the  pressure  increases, 
the  tail  seal  at  the  end  of  the  projectile  flexes,  allowing  some  liquid 
Injection.  At  the  desired  shot  start  pressure,  the  shot  start  link 
breaks,  releasing  the  projectile.  The  projectile  then  rapidly  opens  the 
injection  holes,  leading  to  rapid  propellant  injection.  Near  the  end  of 
the  piston  stroke,  the  piston  partially  blocks  the  injection  ports.  The 
liquid  pressure  rises,  and  brings  the  piston  to  a  gradual  stop. 

Ill .  BASIC  ASSUMPTIONS 

The  propellant  in  the  liquid  reservoir  is  assumed  to  be  a 
homogeneous,  isothermal  fluid.  Because  of  the  high  pressures  in  a  gun, 
the  liquid  is  considered  to  be  compressible. 

The  piston  and  gun  tube  are  assumed  to  be  infinitely  thin,  that  is, 
liquid  moves  instantaneously  from  the  reservoir  to  the  chamber/ tube . 

The  actual  flow  in  the  orifices  is  ignored.  The  mass  flux  out  of  the 
reservoir  is  normally  computed  assuming  steady  state  Bernoulli  flow. 
Comparisons  with  transient  lumped  parameter,  one - dimens inal ,  and  two- 
dimensional  models  of  the  liquid  reservoir  indicate  that  this  is  a 
reasonable  assumption.  As  a  first  approximation,  the  liquid  is  assumed 
to  combust  instantaneously  as  soon  as  it  enters  the  chamber  or  the  gun 
tube,  releasing  all  its  energy.  Later  finite  burning  rate  models  for 
the  liquid  propellant  will  be  introduced,  in  an  attempt  to  model  at 
least  crudely  the  effects  of  propellant  accumulation. 

The  intermediate  chamber  fluid  is  assumed  to  be  homogeneous  and 
stagnant.  Lacking  any  detailed  information  on  flow  in  the  chamber,  the 
fluid  velocity  in  the  chamber  is  simply  set  uniformly  equal  to  zero. 

The  primer  is  assumed  to  have  the  same  thermochemical  properties  as 
the  propellant.  The  combustion  and  injection  of  the  primer  need  not  be 
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considered  in  detail.  Instead  the  initial  pressure  in  the  charaber/tube 
is  an  input  parameter.  Then  the  amount  of  propellant  needed  to  produce 
this  pressure  is  calculated,  and  this  is  taken  as  the  mass  of  the 
primer.  The  piston  and  projectile  are  not  allowed  to  move  during  the 
burning  of  the  primer.  Later  more  complicated  primer  options  will  be 
discussed. 

Unlike  solid  propellant  guns,  there  is  usually  a  large  area  change 
between  the  intermediate  chamber  and  the  gun  tube.  The  mass  flux  into 
the  gun  tube  is  computed  assuming  steady  state  isentropic  flow  without 
loss  terms. 

The  gun  tube  is  also  treated  as  a  lumped  parameter  region,  using 
the  Lagrange  pressure  distribution  that  has  successfully  modeled  the 
behavior  of  solid  propellant  guns® .  The  flow  between  the  chamber  and  the 
tube  is  ignored,  and  the  velocity  at  the  end  of  the  gun  tube  is  set 
equal  to  zero.  Because  of  the  large  area  between  the  chamber  and  the 
tube,  the  flow  velocities  will  be  small  here.  The  gun  tube  is  treated 
as  a  simple  cylinder,  and  the  structure  at  the  breech  end  of  the  gun  is 
ignored.  The  projectile  is  also  treated  as  a  cylinder,  and  the  tail 
seal  is  not  modeled. 

Below,  we  derive  the  governing  equations  for  the  simple  lumped 
parameter  model  described  here.  Then  a  list  of  options  will  be  given, 
with  the  appropriate  modifications  to  the  base  governing  equations. 

IV.  GOVERNING  EQUATIONS 

With  the  above  assumptions ,  the  regenerative  gun  behavior  can  be 
modeled  by  13  ordinary  differential  equations.  These  equations  are 
marked  using  numbers  in  square  brackets  at  the  left  ([1]  to  [13]).  Some 
algebraic  equations  are  required  to  compute  the  coefficients  of  the 
ordinary  differential  equations .  The  equations  describe  the  three 
lumped  parameter  regions  (liquid  reservoir,  intermediate  chamber,  and 
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gun  tube),  the  mass  flux  between  the  regions,  and  the  piston  and 
projectile  motion.  For  historical  reasons,  the  regions  are  numbered  1, 
2,  and  4.  In  a  previous  model  the  back  end  of  the  gun  tube  was  called 
the  combustion  chamber  (region  3).  Since  this  is  just  an  extension  of 
the  gun  tube,  it  is  not  considered  a  separate  region  in  this  report. 
Most  of  the  derivations  are  the  same  as  for  the  in  line  regenerative 
gun,  and  are  not  repeated  in  detail  here.  The  notation  is  given  in  the 
glossary.  For  this  section  the  chamber  and  the  tube  are  assumed  to 
contain  only  gas  phase  products. 


A.  LIQUID  RESERVOIR 


The  equations  governing  the  piston  motion  are  straightforward.  Let 

be  the  volume  of  the  propellant  chamber,  and  A^  be  the  area  of  the 

propellant  side  of  the  piston  (this  area  includes  the  vent  holes).  Let 

s__  be  the  piston  travel  (to  the  right)  and  v__  be  the  piston 
ps  ps 

velocity.  Then 


dt 


-v 


ps 


A1 


(1) 


The  acceleration  of  the  piston  equals  the  force  on  the  piston 
divided  by  the  mass  of  the  piston.  That  is, 


dv  g 

(2]  — —  “  -  [P2  (A2'Ag'A12^  *  (Pl+Pps^  ^A1*A12^  ’  (2) 

dt  MpS 

where  A2  is  the  area  of  the  combustion  side  of  the  piston,  Ag  is  the 
area  of  the  grease  dyke  between  the  piston  and  the  chamber  wall,  A^  is 
the  area  of  the  hole  in  the  piston,  p^  is  the  pressure  in  the  liquid 
chamber,  p£  is  the  pressure  in  the  combustion  chamber,  and  M^s  is  the 

mass  of  the  piston.  The  quantity  gQ  -  10^  gm/s-cm-MPa  is  a  conversion 

o 

constant  to  put  the  acceleration  in  the  desired  units  of  cm  /s.  By 
assumption,  the  pressures  are  constant  throughout  the  two  chambers.  The 
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pressure  in  the  grease  dyke  is  normally  very  close  to  the  chamber 
pressure,  so  it  balances  the  chamber  pressure  on  this  part  of  the 
piston.  The  quantity  p__  is  an  empirical  correction  to  simulate  the 

pa 

effects  of  frictional  resistance  (see  options).  The  piston  is  assumed 
initially  to  be  prevented  from  moving  toward  the  breech.  So  if  the 
initial  acceleration  is  negative,  it  is  changed  to  zero.  Also,  the 
quantity  Ppg  is  solely  a  resistance  pressure.  If  the  Ppg  term  is 
initially  larger  than  the  other  pressure  terms ,  the  acceleration  is  set 
equal  to  zero.  Also,  if  at  a  later  time  the  piston  reverses  direction, 
the  sign  of  pps  is  changed.  The  equation  governing  the  piston  travel  is 


ds 


[3] 


ps 


-  V, 


dt 


ps  • 


From  conservations  of  mass , 


(3) 


[4] 


d Pi  Pi  dVi  m^2 

dt  dt  V1  Vx 


W 


where  p-^  is  the  density  of  the  liquid,  m^  is  the  mass  flux  from  the 
reservoir  into  the  chamber,  and  m^^  is  the  mass  flux  from  the  reservoir 
into  the  gun  tube.  The  energy  equation  can  be  expressed  as 


o 

dpx  c,'-  dpi 

[5]  - ,  (5) 

dt  gQ  dt 

where  c^  is  the  speed  of  sound  in  the  liquid. 

Heat  transfer  and  temperature  changes  in  the  liquid  propellant  are 
expected  to  be  small,  and  are  ignored.  Then  the  speed  of  sound  is  given 
by 


C1  -  1/  go  K  /  Pi  , 


(6) 
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where  K  is  the  adiabatic  bulk  modulus.  For  the  common  liquid 
propellants,  this  can  be  fit  very  accurately  by^ 

K  -  +  K2  px  .  (7) 

The  corresponding  equation  of  state  is 

K1  k2 

Pi -  [ (p\/P0)  «  1]  •  (8) 

k2 

The  ordinary  differential  equation  [5]  for  the  pressure  is  not  actually 
required.  Once  the  density  is  known,  the  pressure  can  be  computed  from 
the  equation  of  state  (8). 

Recently  an  equation  of  state  has  been  derived  that  also  includes 

O 

the  temperature  dependence.  This  could  be  used  instead  to  include 
temperature  changes  in  the  liquid.  The  new  equation  predicts  almost  the 
same  pressure  dependence  as  eq.  (8). 

The  mass  flux  from  the  reservoir  into  the  chamber  is  assumed  to  be 
steady  state  Bernoulli  flow 


m12  ”  CD  A12  J  2S0  P\  (Pi  *  P2^  •  (9) 

The  discharge  coefficient  is  an  empirical  correction  to  the  equations 
to  take  into  account  frictional  losses.  Since  the  mass  flux  for  a  gun 
firing  is  not  steady  state,  Cp  may  also  be  used  to  approximate  the  time 
delay  in  reaching  steady  state. 

Similarly,  the  mass  flux  from  the  reservoir  into  the  gun  tube  is 
given  by 


“14  "  CD  a14  J  2 *1  <Pl  *  Ph>  ■ 


(10) 
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where  p^  is  Che  pressure  in  Che  gun  Cube  aC  Che  injeccion  holes  (see 
SecCion  C  below) .  The  same  discharge  coeff icient  is  used  for  boch 
injeccion  processes. 


The  vent  area  must  be  calculated.  A  "Y"  shaped  vent  is. 

assumed  (see  Figure  2).  The  vent  openings  are  considered  to  be 
circular.  First  the  position  of  the  piston  is  calculated.  If  the 
piston  has  reached  Che  injeccion  holes,  it  partially  blocks  the  holes. 
The  area  not  blocked  by  the  piston  is  calculated.  A  similar  procedure 
is  carried  out  on  the  gun  Cube  side.  The  initial  position  of  Che 
projecCile  is  used  to  calculate  how  much  of  Che  piston  side  holes  are 
blocked.  The  vent  area  between  the  reservoir  and  the  Cube  is  taken  to 
be  the  smaller  of  these  two  areas.  If  a  simple  cylindrical  vent  is  used 
instead  of  a  "Y"  vent,  the  area  of  the  second  set  of  holes  in  the  gun 
tube  is  set  equal  to  zero. 

B.  COMBUSTION  CHAMBER 

Now  consider  the  combustion  chamber.  By  analogy  with  eq.  [1], 
dV2 

[6]  - v  A2  ,  (11) 

dt  * 

where  V2  is  the  volume  of  the  combustion  chamber  and  A2  is  the  area  of 
the  combustion  side  of  the  piston  (including  holes).  Similarly,  from 
conservation  of  mass 


dp2  P2  dV2  m^2  -  m2^ 

dt  V2  dt  V2 


(12) 


where  p2  is  the  density  of  the  gas  in  the  combustion  chamber  and  m2^  is 
the  mass  flux  into  the  gun  tube.  The  energy  equation  is 
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[8] 


(13) 


o 

dp2  C2  dp 2  ni^2  (^Ll  *  ^ 

dt  gQ  dt  V2  -  b  M2 

m24  '  hG2^  ^  ^ 

V2  •  b  M2 

where  c2  Is  the  speed  of  sound  in  the  chamber,  hLl  is  the  enthalpy  of 
the  liquid,  h^  Is  the  enthalpy  of  the  chamber  gas,  is  the- enthalpy 
of  the  gun  tube  gas,  and  m2^  is  the  mass  flux  from  the  chamber  into  the 
tube.  The  last  term  only  appears  if  m2^  is  negative,  since  flux  out  of 
the  region  will  not  effect  the  enthalpy.  The  energy  equation  is  based 
on  the  Noble-Abel  equation  of  state® 

P2  -  P2  Rs  t2  /  <1  -  b  *  (14) 

where  T2  is  the  temperature  in  region  2,  Ra  is  the  specific  gas  constant 
(universal  gas  constant  divided  by  molecular  weight  of  the  gas) ,  and  b 
is  the  covolume.  This  can  also  be  written  as 

P2  -  P2  -  1)  cv  T2  /  (1  -  b  p2)  ,  (15) 

where  cv  is  the  specific  heat  at  constant  volume  and  7  is  the  specific 
heat  ratio  Cp/cv.  The  corresponding  equation  for  the  speed  of  sound  is 


So  y  P2 

c  2  -  ■/  - 

p2  (1  -  b  p2) 


(16) 


The  energy  content  of  a  propellant  is  normally  given  in  terms  of 
impetus  A.  Consider  a  quantity  of  propellant  in  a  constant  volume 
closed  chamber.  The  propellant  combusts,  and  the  result  is  a  gas  with 
some  molecular  weight  Mg  and  temperature  TQ  (called  the  isochoric 
temperature) .  The  impetus  is  defined  as 


A  “  *u  To  /  Mg  ’ 


(17) 
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where  ^  is  Che  universal  gas  constant.  Impetus  can  be  obtained  from  a 
closed  bomb  experiment,  if  the  proper  corrections  are  made  for  the 
primer  and  the  heat  loss  to  the  chamber  walls.  For  liquid  propellants, 

Q 

the  impetus  has  been  calculated  by  a  thermodynamics  code.  The  values 
of  the  molecular  weight,  the  specific  heats,  and  the  covolume  of  the 
product  gases  are  evaluated  at  the  isochoric  temperature  and  then 
assumed  to  be  constant  with  respect  to  temperature  and  pressure.  The 
impetus  does  depend  on  the  initial  loading  density  (grams  of  propellant 
per  unit  volume).  Values  are  chosen  at  a  nominal  loading  density  of  0.2 
g/cc.  If  the  values  are  instead  chosen  for  some  other  loading  density, 
the  effects  are  negligable. 

For  our  case,  it  is  more  convenient  to  work  with  the  chemical 
energy 


ei  -  A  /  (7-I)  .  (18) 

The  enthalpy  of  the  liquid  is  given  by 

hLl  "  ®l  +  Pi  /  •  (19) 

The  enthalpy  of  the  gas  is  given  by 

^G2  "  cv  ^2  +  P2  /  p2  “  cp  ^2  +  k  P2  •  (20) 

C.  GUN  TUBE 


Finally,  the  gun  tube  (region  4)  is  considered  to  be  a  simple 
cylinder.  The  structure  shown  in  Figure  2  is  ignored.  The  volume  has 
the  standard  equation 
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[9] 


Vpj  A4  ’ 


(21) 


dV, 


dc 


where  Vpj  is  the  velocity  of  the  projectile  and  is  the  area  of  the 
gun  tube.  The  rapid  projectile  motion  creates  a  large  pressure 
gradient.  The  standard  approach  is  to  use  a  Lagrange  pressure 
distribution, ^  that  is,  assume  that  the  density  is  constant  with  respect 
to  space.  It  follows  from  the  one -dimensional  continuity  equation  that 
the  gas  velocity  is  a  linear  function  of  distance.  The  gas  velocity  at 
the  base  of  the  projectile  is  equal  to  the  velocity  of  the  projectile. 

In  solid  propellant  guns,  the  gas  velocity  at  the  other  end  of  the  gun 
(breech)  is  assumed  to  be  zero.  Since  there  is  a  mass  flow  into  the  gun 
tube,  this  is  not  correct  for  our  model.  However,  present  designs  have 
a  large  area  between  the  chamber  and  the  tube,  and  the  corresponding 
velocity  is  relatively  small.  In  this  case  the  standard  Lagrange 
aproximation  is  accurate. 


Integrating  the  momentum  equation, 


M4  x2 

P (x)  -  PL  *  (Pr  *  Ppj) -  — j  •  (22> 

2Mpj  XR 

where  xR  is  the  distance  from  the  back  end  of  the  tube  to  the  base  of 
the  projectile,  p^  is  the  pressure  at  the  back  of  the  gun  tube,  pR  is 
the  pressure  at  the  base  of  the  projectile,  is  the  mass  of  the  gas  in 

the  tube,  Mpj  is  the  mass  of  the  projectile,  and  Ppj  is  the  resistance 
pressure.  The  latter  is  an  input  parameter,  which  takes  into  account 
the  shot  start  engraving  force  and  the  frictional  forces  between  the 
projectile  and  the  bore.  Evaluating  at  x-xR  yields 

pR  -  p  -  (p  .  Ppj)  •  (23) 

2MPJ 

Integrating  equation  (21)  from  the  breech  to  the  projectile  base  yields 
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(24) 


M, 


P4  “  PL 


6MPJ 


<PR  *  Ppj>  • 


where  p4  is  the  space  mean  pressure.  Then  the  above  two  equations  can 
be  solved  for  p^  and  p^.  Equation  (22)  can  be  used  to  compute  p^,  which 
is  taken  to  be  the  pressure  at  the  center  of  the  first  set  of  gun  tube 
holes . 


The  above  two  equations  assume  that  the  projectile  has  started  to 
move,  creating  a  pressure  differential.  If  the  projectile  has  not 
moved,  then  region  4  is  treated  as  a  stagnant  homogenous  region  like 
region  2,  and  the  pressure  is  uniform. 

Given  the  pressure  on  the  projectile,  the  acceleration  equation  is 


HQ]  —  -  (pr  -  Ppj)  a4  go  /  Mpj  . 

dt 


(25) 


As  with  the  piston,  the  pressure  Ppj  is  only  a  resistive  pressure.  If 
Ppj  >  PR’  c^e  Pr°jectHe  does  not  move.  The  projectile  travel  Spj  is 
given  by 

d3pj 

[11]  ~ “  -  v  ,  .  (26) 

dt 

The  last  two  differential  equations  are  exactly  analogous  to  region  2 . 
That  is, 


dp4  Pl.  <1^4  “14  “94 

[12]  - +  -  +  -  .  (27) 

dt  v4  dt  v4  v4 


where  p4  is  the  space  mean  average  density  of  the  gas  in  region  4,  and 
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(28) 


[13]  ^  °4  dp4  +  mi4  ^  '  h°4>  ^  ’  1) 

dt  gQ  dt  V4  -  b  M4 


“24  ^2  *  ^34^  ^ 


V4  -  b  M4 

The  last  term  only  appears  If  m24  is  positive  (p2  >  p^) •  The  average 
speed  of  sound  in  the  tube  is 


~  J 


So  T  P4 

P4  (1  -  b  p4) 


The  mass  flux  from  the  chamber  into  the  gun  tube  is  assumed  to  be 
isentropic  flow  (that  is,  adiabatic  and  reversible) . ^  Assume  that 
p2  is  greater  than  p^.  Then  m24  is  positive.  The  gas  in  the  chamber  is 
considered  to  be  stagnant,  and  expands  isentropically  into  the  gun  tube. 
The  expanded  gas  then  mixes  with  the  gun  tube  gas  already  present.  The 
process  equations  for  a  Noble-Abel  equation  of  state  are 


T  (1/p  -  b)(T*l)  „  constant  , 


p  (1/p  -  b)T  -  constant 


Assume  that  pt  -  p^,  since  the  pressure  will  equilibrate  much  more 
rapidly  than  the  temperature  or  density.  The  process  equations  (30)  and 
(31)  can  be  used  to  find  the  throat  temperature  Tt  and  density  pt.  Then 
the  one -dimensional  momentum  equation  can  be  integrated  from  the 
stagnant  conditions  in  the  chamber  to  the  throat,  resulting  in 


Vt  -  y  2g0  {b(p2-pt)  +  c  Tt  [(P2/Pt)(7~1)/7  -  1]  (32) 


a; 


“24  ”  Cd'  a24  pt  vt  > 


(33) 


where  *-s  c^e  area  between  the  chamber  and  the  tube.  The  discharge 
coefficient  CD'  is  as  before  an  empirical  correction  for  loss  terms. 
Due  to  lack  of  better  information,  this  is  set  to  one. 

If  the  flow  is  from  the  tube  into  the  chamber,  then  "  P2‘ 
process  equations  (30)  and  (31)  are  again  used  to  find  the  throat 
temperature  Tp  and  density  pt.  The  flow  is  described  by 

vc  -  J  2gc  (b(pL-pt)  +  cpTt  [(pL/pt)(7-1)/^  -  1]  (34) 


and 


m24  “  -V  a24  pt  vt  • 

V.  VENT  OPTIONS 

Unlike  the  in-line  regenerative  gun  code,  there  is  a  standard 
setup  for  the  vents.  The  vent  area  between  the  reservoir  and  the 
chamber  is  constant. 

The  vents  into  the  gun  tube  are  "Y"  shaped.  On  the  reservoir 
side,  there  is  a  single  set  of  circular  holes.  The  area  of  each  hole 
and  the  number  of  holes  is  read  in,  as  well  as  the  distance  the  piston 
must  travel  to  reach  the  center  of  the  holes.  As  the  piston  moves,  a 
check  is  made  to  see  if  the  piston  partially  blocks  the  holes.  If  so, 
the  area  covered  by  the  piston  is  computed.  On  the  gun  tube  side, 
there  are  two  sets  of  holes.  The  gun  tube  is  treated  as  a  simple 
cylinder.  The  area  of  each  set  of  holes  is  read  in,  as  well  as  the 
distance  from  the  end  of  the  gun  tube  to  the  center  of  the  holes.  The 
number  of  holes  is  assumed  to  match  the  number  of  holes  on  the  reservoir 
side.  The  projectile  is  also  treated  as  a  simple  cylinder.  The  initial 
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offset  of  the  projectile  from  the  end  of  the  tube  is  entered.  From  the 
position  of  the  holes  and  the  projectile,  the  area  blocked  by  the 
projectile  is  computed.  The  actual  vent  area  A^  is  taken  to  be  the 
minimum  of  the  reservoir  and  tube  areas.  There  is  no  attempt  to  model 
the  tail  seal. 

In  practice,  since  the  holes  in  the  short  barrel  are  drilled  on  a 
slant,  the  openings  will  be  elliptical  rather  than  circular.  Using 
circular  openings  in  the  code  results  in  only  minor  errors  in  the  vent 
area. 

VI.  PISTON  RESISTANCE  OPTIONS 

As  the  piston  moves,  there  is  frictional  resistance  between  the 
piston  and  the  chamber  walls.  There  is  no  good  estimate  for  the  size  of 
this  frictional  resistance.  Because  of  the  large  pressures  involved  and 
the  relatively  slow  'speed  of  the  piston,  this  effect  is  not  expected  to 
be  important.  Nevertheless,  piston  resistance  options  are  included. 

Only  one  option  is  presently  in  the  code. 

A.  PIS1 

The  resistance  is  considered  to  be  solely  a  function  of  piston 
travel.  A  table  of  piston  travel  versus  resistive  pressures  is  read  in. 
The  table  will  be  normalized  to  the  maximum  piston  travel.  The  piston 
resistive  pressure  pps  at  any  given  piston  travel  sps  is  found  by 
interpolation.  The  actual  resistive  force  is  pps  (A^  -  Aj^).  An 
equivalent  pressure  is  read  in  rather  than  the  actual  force  to  make  it 
easier  to  compare  the  braking  force  exerted  by  friction  with  the  braking 
force  exerted  by  the  liquid. 

VII.  PROJECTILE  RESISTANCE  OPTIONS 

Similarly,  as  the  projectile  moves,  there  is  friction  between  the 


15 


projectile  and  the  gun  tube.  More  important  is  the  shot  start  pressure. 
The  projectile  is  restrained  from  moving  until  the  pressure  at  its  base 
reaches  an  appropriate  value.  The  choice  of  shot  start  pressure  can 
have  a  large  effect  on  performance  and  on  the  maximum  pressures  in  the 
gun.  At  present  there  is  only  one  option. 

A.  PR0J1 

As  with  the  piston,  a  table  of  projectile  travel  versus  resistive 
pressure  is  read  in.  The  table  is  not  normalized.  The  resistive 
pressure  Ppj  for  any  given  projectile  travel  Spj  is  found  by 
interpolation.  If  Spj  is  larger  than  the  last  table  entry  for  travel, 
the  last  table  entry  for  pressure  is  chosen.  The  actual  resistive 
force  is  ppj  A^.  The  first  entry  in  the  table  is  the  shot  start 
pressure.  The  projectile  will  not  move  until  the  tube  pressure  reaches 
the  shot  start  pressure. 

VIII.  DISCHARGE  COEFFICIENT  -  FROM  RESERVOIR 

The  discharge  coefficients  needed  in  equations  (9)  and  (10)  are 
input  parameters.  The  same  discharge  coefficient  is  used  for  both 
injection  processes  out  of  the  reservoir. 

A.  DIS1 

A  table  of  piston  travel  versus  discharge  coefficient  is  read  in. 
The  table  is  normalized  to  the  maximum  piston  travel.  The  discharge 
coefficient  Cq  for  any  given  piston  travel  SpS  is  found  by 
interpolation. 

IX.  DISCHARGE  COEFFICIENT  -  FROM  CHAMBER 

The  discharge  coefficient  needed  in  equations  (33)  and  (35)  is 
assumed  to  be  one.  There  is  normally  a  large  area  between  the 
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intermediate  chamber  and  the  gun  tube,  and  loss  terms  are  expected  to  be 
very  small. 


X.  MASS  FLUX  OPTIONS  -  FROM  RESERVOIR 


A.  FLUX1 

Steady  state  Bernoulli  flow  is  assumed.  The  mass  flux  m^  is 
computed  using  equation  (9),  and  m^  is  computed  using  equation  (10). 
The  routine  computes  the  Weber  number,  the  Reynold's  number,  and  the 
third  Lagrange  number  for  m^,  which  is  the  primary  injection  flow. 
While  this  code  does  not  use  these  numbers ,  they  are  printed  out  as 
additional  information  about  the  injection  process. 

B.  FLUX2 

Since  the  gun  environment  is  highly  transient,  steady  state  flow 

may  not  be  a  good  approximation.  A  transient  model  for  injection  has 
5  12 

been  developed.  ’*•  At  the  present  this  is  only  implemented  for  flow 
m^.  For  this  option  the  vent  area  in  the  piston  A^  must  be  set  to 
zero  (no  Injection  through  the  piston  face).  The  basic  assumption  is 
that  the  space  derivative  of  the  mass  flux  pvA  of  the  fluid  in  the 
reservoir  is  zero.  Then  the  one -dimensional  momentum  equations  can  be 
integrated  from  the  piston  to  the  vent  with  the  result 

d(pVA)  22  2  f 

[FI] - [°-5Pl<vps-v14  )  +  g0CD  <Pl'Ph^/Jdx/A  (36) 

dt 

where  v^  is  the  velocity  of  the  fluid  entering  the  tube.  This 
additional  ordinary  differential  equation  is  integrated  to  obtain  the 
mass  flux  rate  into  the  chamber.  The  equation  shows  a  rapid  rise  to 
values  near  steady  state. 
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XI.  MASS  FLUX  OPTIONS  -  FROM  CHAMBER 


A.  FLUX1 . 

Steady  state  isentropic  flow  is  assumed.  The  mass  flux  into  the 
gun  tube  is  computed  using  eq.  (33)  or  (35). 

XII.  GUN  TUBE  PRESSURE  DISTRIBUTION  OPTIONS 

Recently  a  number  of  different  assumptions  for  the  gun  tube 

pressure  distribution  were  implemented  and  tested  against  a  one- 
13 

dimensional  code.  However,  for  the  present  case  the  standard  Lagrange 
distribution  is  reasonable. 

A.  TUBE1 

The  standard  Lagrange  pressure  distribution  is  assumed,  that  is, 
the  density  is  constant  in  the  tube,  which  implies  that  the  velocity  is 
linear.  For  the  purpose  of  computing  the  pressure  distribution,  the 
velocity  at  the  back  end  of  the  gun  tube  is  set  equal  to  zero.  The 
momentum  equation  is  integrated  to  obtain  a  formula  for  the  pressure 
p(x)  in  the  tube.  The  pressure  p^  at  the  base  of  the  projectile  and 
the  pressure  p^  at  the  back  end  of  the  gun  tube  are  computed  using  eq. 
(23)  and  (24).  The  injection  into  the  gun  tube  will  effect  the  pressure 
distribution.  However,  the  injection  velocities  are  small  compared  to 
the  projectile  velocity.  Since  the  injection  is  radial,  there  is  no 
obvious  way  to  modify  the  Lagrange  distribution  to  take  the  injection 
into  account. 


XIII.  PRIMER  OPTIONS 

In  present  liquid  propellant  guns,  a  solid  primer  is  used.  The 
primer  is  in  the  booster  cavity.  The  primer  injects  hot  gas  through 
large  holes  into  the  gun  tube/chamber  to  begin  the  firing  cycle.  The 


18 


exact  details  of  the  primer  model  have  only  a  small  effect  on  pressure, 
so  simplifications  have  been  made.  First,  the  primer  is  assumed  to  be 
the  same  material  as  the  liquid  propellant,  so  a  new  set  of  properties 
for  the  actual  solid  primer  does  not  have  to  be  read  in.  Three  primer 
options  then  offer  varying  levels  of  simplification. 

A.  PRIM1 

It  is  assumed  that  the  primer  completely  combusts  before  the  start 
of  the  calculation,  pressurizing  the  chamber.  The  corresponding  piston 
movement  is  ignored.  The  code  computes  how  much  propellant  would  be 
necessary  to  create  the  input  chamber  pressure,  and  records  this  as  the 
primer  weight. 

B.  PRIM2 

The  initial  mass  of  the  primer  is  read  in.  Then  the  primer  is 
distributed  evenly  in  the  chamber  and  tube  as  liquid  droplets.  A 
droplet  burning  option  must  be  chosen  (see  below).  The  size  of  the 
droplets  is  determined  from  the  drop  diameter  in  the  droplet  burning 
option.  This  option  can  be  used  to  mimic  the  delay  in  time  to  reach  the 
primer  pressure. 

C.  PRIM3 

The  primer  is  out  of  the  system.  The  mass  of  the  primer,  the 
injection  time,  and  a  heat  loss  term  are  read  in.  The  primer  is  then 
injected  at  a  steady  rate  into  the  gun  tube  over  the  desired  time 
interval.  The  primer  is  assumed  to  combust  instantaneously  as  it 
enters  the  tube.  Since  primer  gases  are  often  injected  through  long 
narrow  tubes,  heat  loss  is  important.  To  take  heat  loss  into  account, 
the  chemical  energy  of  the  propellant  is  multiplied  by  the  heat  loss 
factor.  During  the  injection  process,  the  gun  tube  density  and  pressure 
equations  are  modified 
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[12] 


(37) 


d/>4 

dC 


[13] 


dp4 

dt 


V  (VS  -  (7*D 

+  - - - -  . 

V4  -  b  H4 


(38) 


where  nip4  Is  the  mass  flux  of  primer  into  the  tube  and  lp  is  the  heat 
loss  term. 


XIV.  GUN  TUBE  HEAT  LOSS  OPTIONS 


In  gun  systems  there  is  a  large  energy  flux  from  the  hot  gas  to 
the  gun  tube.  This  is  important  for  actual  systems  primarily  because  it 
causes  erosion  of  the  gun  tube.  The  effect  on  muzzle  velocity  and 
pressure  is  on  the  order  of  a  few  percent.  Heat  loss  is  often  treated 
as  an  adjustable  parameter  for  fine  tuning  a  model  to  agree  with 
experiment. 

A.  HEAT1 


The  heat  loss  to  the  gun  tube  walls  is  ignored.  This  causes  a 
slight  increase  In  muzzle  velocity  and  maximum  pressures. 

B.  HEAT2 

There  is  heat  loss  to  the  gun  tube.  A  correlation  developed  for 

O 

turbulent  pipe  flow  is  used.  An  assumed  gun  tube  temperature  Tw  is 
read  in  and  assumed  to  remain  constant.  The  equation  for  the  space  mean 
pressure  is  modified  by 
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(131 


(39) 


dp4  qw  (7  -  1) 

dt  (1  -  b  pu) 

where  Qw  is  Che  heat  loss  Co  Che  Cube  wall.  SimulaCions  run  so  far  wich 
Che  code  indicaCe  a  heac  loss  of  3%  Co  5%  for  120mm  guns,  and  10%  Co  12% 
for  30mm  guns.  These  numbers  are  in  Che  range  reporced  by  Nordheim, 
Soodak,  and  Nordheim^  for  Cheir  analysis  of  experimenCal  daca,  and 
comparable  wich  Che  heaC  loss  normally  used  in  solid  propellanc  gun 
modeling.  If  desired,  Che  heac  loss  can  be  adjusced  up  or  down. 

XV.  AIR  SHOCK 

The  projeccile  will  compress  Che  air  in  Che  Cube  as  ic 
acceleraCes.  This  will  have  a  minor  effect  on  Che  velocicy  of  Che 
proj eccile. 

A.  SH0CK1 

The  resisCance  due  Co  Che  air  shock  in  fronc  of  Che  projeccile  is 
ignored. 

B.  SH0CK2 

The  resisCance  pressure  p_  is  calculaCed  using  a  formula  developed 
by  Corner**  for  Che  air  shock  based  on  Che  Rankine-HugonioC  relacions. 
This  new  resisCance  pressure  ps  is  jusc  added  on  Co  Che  fricCional 
resisCance  pressure  Ppj . 


XVI .  REPEAT  OPTIONS 

In  sysCems  scudies,  Chere  is  ofCen  a  physical  conscrainC  on  Che 
class  of  accepCable  solucions.  One  common  conscrainC  is  implemenced  in 
che  code.  The  code  will  inCegraCe  repeaCedly  unCil  Che  conscrainC  is 
saCisf led. 
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A.  REP1 


The  code  integrates  once.  This  is  the  standard  option. 

B.  REP2 

A  solution  is  sought  with  a  given  maximum  liquid  pressure.  The 
liquid  pressure  is  adjusted  by  changing  the  vent  area  into  the  gun  tube 
on  the  liquid  side.  The  desired  liquid  pressure  and  a  first  guess  for 
the  vent  increment  is  read  in. 

The  problem  is  first  integrated  with  the  input  conditions.  If  the 
maximum  liquid  pressure  is  too  low,  the  vent  area  is  incremented  by  the 
input  value.  If  the  liquid  pressure  is  too  low,  the  vent  area  is 
decremented.  Once  values  are  obtained  both  above  and  below  the  desired 
maximum  liquid  pressure,  interpolation  is  used  to  obtain  the  new  vent 
area.  Convergence  normally  occurs  after  a  small  number  of  integrations. 

XVII.  CHAMBER  PRESSURE  OPTIONS 

It  is  convenient  to  isolate  part  of  the  model  for  study.  For 
instance,  suppose  that  the  model  chamber  pressure  Is  slightly  high. 

This  raises  the  piston  acceleration  and  increases  the  rate  of  injection 
of  the  propellant.  As  the  propellant  bums,  the  chamber  pressure  rises. 
After  a  short  time,  the  model  chamber  pressure  and  piston  travel  are 
substantially  too  large.  So  to  study  the  liquid  injection  process,  it 
is  better  to  disconnect  the  chamber  pressure  from  the  injection. 

A.  CHAM1 

The  chamber  pressure  is  computed  using  the  usual  equations. 
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B.  CHAM2 


The  experimental  chamber  pressure  versus  time  is  read  in.  At  each 
time  step,  an  interpolated  experimental  chamber  pressure  overwrites  the 
pressure  computed  using  the  ordinary  differential  equation  [8].  This 
allows  a  study  of  the  injection  process  and  the  projectile  motion  using 
the  experimental  data  as  a  boundary  condition. 

C.  CHAM3 

The  time  derivative  of  the  experimental  chamber  pressure  is  read 
in.  This  is  normally  found  by  making  a  spline  fit  of  the  chamber 
pressure  and  taking  the  derivative  of  the  spline.  At  each  time  step,  an 
interpolated  derivative  overwrites  the  derivative  computed  from  the 
ordinary  differential  equation  [8].  This  is  more  stable  than  the  CHAM2 
option,  and  larger  time  steps  will  be  taken.  The  code  runs 
substantially  faster,  although  the  same  answer  will  be  obtained. 

XVIII.  DROPLET  BURNING  OPTIONS 

So  far  the  liquid  has  been  assumed  to  combust  instantaneously  upon 
exiting  the  reservoir.  However,  liquid  accumulation  in  the  gun  tube  is 
very  important  for  some  cases.  So  simple  rules  are  derived  for  the 
formation  and  combustion  of  droplets  in  the  gun  tube/chamber.  In  the 
model,  the  liquid  is  assumed  to  instantaneously  form  droplets  of  fixed 
size  as  the  liquid  enters  the  gun  tube.  The  droplets  are  evenly 
distributed  in  the  tube.  Flow  into  the  intermediate  chamber  will 
include  droplets.  The  droplets  will  then  combust  at  some  given  rate. 

The  initial  size  of  the  droplets  is  an  input  parameter.  This  simplified 
model  allows  the  code  to  determine  the  effects  of  liquid  accumulation. 

The  derivation  of  the  relevant  equations  is  very  similar  to  the 
previous  reports.  Let  Mj^  be  the  liquid  mass  in  region  2-  VL2  be  the 
liquid  volume,  Mq2  be  the  gas  mass,  and  VG2  be  the  gas  volume .  Then 
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(40) 


*L2 

PL2 - 

VL2 


and 


*02 

Pq2  *  .  (41) 

VL2 


where  p ^  is  the  density  of  the  liquid  and  p G2  is  the  density  of  the 
gas.  Define  the  porosity 


Then 


(42) 


P2 


*e2>  PL2  +  e2  ^G2 


(43) 


Following  eq.  (6),  the  speed  of  sound  in  the  liquid  is 

CL2  "  J  So  K  /  PL2  ’  (44) 

and  following  eq.  (16)  the  speed  of  sound  in  the  gas  is 


.  So  ^  P2 

cG2  -  J - 

Pq2  (1  -  b  Pq2) 


(45) 


Assuming  the  gas  and  liquid  pressure  are  the  same,  the  speed  of  sound  c2 
in  the  mixture  is  given  by 
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1 


(46) 


1 

c2  “  J  1  2  2 

P2  «2/^pG2  CG2  )  +  C1  *  «2)/^L2  cL2  ) 

The  enthalpy  of  the  liquid  is 

hL2  “  el  +  p2  /  PH  7 ) 

where  e^  is  the  chemical  energy  of  the  liquid,  and  the  enthalpy  of  the 
gas  is 

hG2  “  cp  T2  +  b  p2  •  (48) 


The  thirteen  equations  derived  in  Section  IV  remain  the  same 
except  for  the  energy  equations  [8]  and  [13].  The  energy  equation  is 
written  in  the  more  general  form 


[8] 


2 

dp2  f>2  c2 
So  V2 


dV2  ^  mL2  ^  nG2  | 


dt 


p L2  pG2 


(49) 


S2  So  *  1) 

- 2 -  }  * 

pG2  cG2  (^  *  b  PQ2 ) 

Here  m^2  is  the  rate  of  production  of  liquid  in  the  chamber.  In  this 
case , 

m24>0  (50) 

m24<0  (51) 


M- 


m^2  “  ni^2  *  ~ 


•L2 
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M. 


“L2  “  m12  *  “2  * 


L4 
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M, 


with  m2  the  rate  at  which  combustion  is  turning  the  liquid  into  gas. 
The  rate  depends  on  whether  fluid  is  flowing  from  the  chamber  to  the 


tube  (m24>0)  or  from  the  tube  into  the  chamber  (m2^<0).  The 
production  of  gas  m^  is  given  by 


“G2 


* 


m24>0 


®G2 


®24<® 


The  rate  of  production  of  mass  times  the  change  in  enthalpy 
by 


S2  ~  m2  (h^2  -  hg2)  +  m^2  (h^  -  h^2)  ,  m2^>0 

^2  "  ®2  ^bL2  "  bG2^  +  m12  (bLl  *  bL2^  ®24<“® 

^4  ML4 

-  m24  (hQ4  -  hG2)  -  -  m24  (h^  -  h^)  -  . 

m4  ■  m4 

In  the  gun  tube,  the  equations  are  exactly  analogous, 
difference  is  in  the  coefficients  for  the  pressure  equation 


2 

dP4  P 4  C4  m^  m^ 

[13] - { - +  -  +  -  + 

dt  S0  v4  dt  ?L4  ^G4 

s4  S0  (7  ■  1) 

- 5 -  )  • 

^G4  CG4  (i  '  b  Pg4^ 

In  the  gun  tube 


mL4  ”  ®14  '  ®4  +  m24  -  >  m24>® 

m2 


rate  of 

(52) 

(53) 

S2  is  given 

(54) 

(55) 

The  only 

(56) 


(57) 
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mL4  “  m14  *  m4  +  m24  — 

m4 


*G2 

“G4  “  m4  +  m24  -  +  °p4  » 

m2 


*G4 

mG4  “  m4  +  m24  -  +  ®p4  ’ 

m4 


m24<®  (58) 


m24>0  (59) 


m24<0  (60) 


S4  “  m4  (h^  "  hg4^  +  m14  ^L1  *  ^L4^  m24'>®  (81) 

^2  ML2 

+  m24  (hG2  -  hG4)  -  +  m24  (hL2  -  hj^) 

M2  M2 

+  v  <ei  •  • 

S4  “  m4  (hj^  -  hG4)  +  m^4  (h^  -  hj^)  m24<0  (62) 

+  “^4  (el  *  hG4^  ’ 

If  there  is  heat  loss  in  the  gun  tube,  the  modification  term  from  eq. 
(39)  is  included. 

To  complete  the  system,  ordinary  differential  equations  are 
required  to  determine  the  liquid  density  and  mass.  The  corresponding 
gas  quantities  can  then  be  easily  derived.  The  basic  assumption  is  that 
the  pressure  in  the  liquid  and  the  gas  in  a  given  region  is  the  same. 
Since  the  gas  pressure  is  known,  the  liquid  equation  of  state  can  be 
written  in  differential  form  as 


dpL2  80  d?2 

[14] - -  - 

dt  cL2  dt 


(63) 
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and 


[16] 


dfilA  S0  dP4 
dt  c dt 


(64) 


The  mass  conservation  equations  for  the  liquid  have  already  been  derived 


[15] 


mL2  ■ 


and 


(65) 


[17] - .  (66) 

dt 

To  close  the  system,  the  rate  at  which  the  liquid  droplets  are 
combusting  must  be  known.  The  rate  of  surface  regression  is  assumed  to 
be  of  the  form 


rate  of  surface  regression 


(67) 


Liquid  propellant  burning  rates  have  been  measured  by  McBratney^ ’ ^ 
The  rate  of  combustion  is 


m2  -  pL2  S  Ap®  ,  (68) 

where  S  is  the  total  surface  area  of  the  droplets  in  a  region.  All  the 
droplets  are  assumed  to  have  a  constant  diameter  d.  So  the  total 
surface  area  in  the  region 


S  -  6  VL2  /d  . 

and  eq.  (68)  can  be  written  as 


(69) 
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(70) 


m2  “  ML2  <6/d)  a?2B  • 

For  region  4,  the  pressure  is  no  longer  constant  over  the  region, 
but  follows  the  Lagrange  pressure  distribution.  This  would  be  quite 
complicated  to  keep  track  of,  so  the  combustion  rate  is  assumed  to 
depend  on  the  average  pressure  p4  ,  and 

o4  -  Mu  (6/d)  Ap4B  .  (71) 

The  droplet  options  can  now  be  described. 

A.  DR0P1 

The  liquid  is  assumed  to  combust  instantaneously  as  soon  as  it 
enters  the  gun  tube/chamber.  This  is  the  model  described  in  section  IV. 

B.  DR0P2 

The  liquid  exits  the  reservoir  and  instantaneously  forms  droplets 
of  diameter  d.  The  size  of  the  droplets  is  fixed  in  time  and  space.  As 
the  droplets  bum,  the  number  of  droplets  may  decrease,  but  the  size  of 
the  droplets  remains  unchanged. 

The  droplet  diameter  and  the  burning. rate  coefficients  are  read 
in.  There  is  some  evidence  that  the  burning  rate  for  liquid 
monopropellants  has  a  sharp  slope  break  around  100  MPa.  So  the  burning 
rate  is  read  in  as  two  functions,  and  the  pressure  is  specified  at  which 
the  change  is  made  from  the  low  pressure  rate  to  the  high  pressure  rate . 

C .  DR0P3 

The  droplet  size  is  fixed  in  space,  but  not  in  time.  A  table  of 
droplet  diameters  versus  piston  travel  is  read  in.  As  usual,  the  table 
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is  scaled  Co  Che  maximum  piscon  Cravel. 

While  Chis  model  is  unrealiscic,  ic  allows  Che  modeling  of  various 
accumulacion  races  ac  differenC  cimes  in  Che  firing  cycle. 


D.  DR0P4 

This  is  like  Che  opcion  DR0P3,  excepC  Che  Cable  of  droplec 
diameCers  is  read  in  versus  Che  projeccile  Cravel. 

XIX.  MASS  AND  ENERGY  BALANCE 

As  a  check  on  Che  equacions  derived,  boch  mass  and  energy  should 
be  conserved.  The  mass  balance  is  scraighcforward. 

MT  -  Mj  +  +  Mq2  +  +  Mp  ,  (72) 

where  Mp  is  Che  mass  of  Che  primer  and  Mj  is  Che  cocal  mass  of  Che 
propellanc  and  primer.  M^.  should  be  consCanC  chroughouc  Che 
incegracion. 

The  energy  balance  is  more  complicaCed.  The  liquid  is  considered 
Co  be  an  isochermal  fluid,  and  iCs  inCemal  energy  e^  is  jusc  Che 
chemical  energy  of  Che  propellanc.  The  CoCal  energy  in  region  1  is 

Ex  -  e1  Mj_  .  (73) 

If  che  oucside  primer  opcion  is  chosen  (PRIM3)  Chen 

Ep-eiMp.  (74) 

The  incernal  energy  of  che  liquid  in  region  2  is  given  by 
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EL2  -  el  ML2  ’  (75> 

and  the  internal  gas  energy  is 

EG2  "  Cv  T2  ^G2  •  <76) 

Similarly  for  region  4 

eL4  -  el  mL4  •  (77> 

and 

EG4  “  cv  T4  «G4  •  (78) 

The  kinetic  energy  of  the  piston  is  given  by 

BKp,  -  °’5  V  vps2  '  So  '  <79> 

and  the  kinetic  energy  of  the  projectile  by 

SKpj  -  0.5  Mpj  vpj2  /  go  .  (80) 

The  liquid  and  gas  in  region  4  are  moving,  and  are  assumed  to  have 
the  same  velocity.  The  kinetic  energy  of  the  liquid  is 

EKl A  ’  °-5  *L4  vpj2  XR  /(3go>  -  0*1) 

where  xR  is  the  distance  from  the  back  of  the  gun  tube  to  the  base  of 
the  projectile.  The  kinetic  energy  of  the  gas  is 

EKG4  -  0.5  Mg4  vpj2  xR  /(3g0)  .  (82) 

The  fluid  in  the  chamber  is  considered  stagnant.  When  the  fluid 
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enters  the  gun  tube,  it  instantaneously  acquires  kinetic  energy.  To 
keep  the  energy  balanced,  the  enthalpy  of  the  liquid  in  the  gun  tube  is 
redefined  as 


hj^  “  +  P4  /  Pl4  +  EKlA  f  **4  *  (83) 

and  the  enthalpy  of  the  gas  in  the  gun  tube  is 

hG4  “  cp  T4  +  b  p4  +  EK^  /  M4  .  (84) 

The  heat  loss  to  the  gun  tube  is  energy  that  leaves  the  system. 

The  total  energy  loss  is  approximated  by  finite  differences.  The  code 
is  set  up  to  print  out  information  at  designated  time  intervals.  After 
each  time  interval,  the  heat  loss  from  the  gun  tube  is  approximated  by 

EH4(t)  -  EH4(t0)  +  Qw  V4  <t-t0)  ,  (85) 

where  Qw  and  V4  are  averages  of  the  values  at  the  present  time  t  and  the 
old  time  tQ. 

The  energy  lost  through  friction  is  approximated  in  a  similar 
fashion.  The  energy  lost  is  equal  to  force  times  distance,  or 

EFps(t)  -  EFps(t0)  +  Pps  (ArA12)(sps(t)-sps(t0))  (86) 

for  the  piston  and 

EFpj (t)  -  EFpJ(t0)  +  Ppj  A4  (spj(t)-spj(t0))  (87) 
for  the  projectile. 

Then  the  total  energy  of  the  system  is 
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®T  "  EL1  +  Ep  +  EL2  +  EG2  +  EL4  +  EG4  +  +  (88> 

EKpj  +  EKu  +  EK^  +  EH4  +  EFps  +  EFpJ  . 

This  should  be  constant  throughout  the  integration.  In  practice,  energy 
errors  of  up  to  0.5%  have  been  observed.  There  is  a  small  error  in  the 
liquid  representation.  As  the  liquid  is  compressed,  the  temperature 
will  rise  slightly,  leading  to  a  larger  internal  energy.  This  effect  is 
ignored  in  the  model.  Other  small  errors  are  introduced  by  the 
approximations  used  for  the  heat  and  frictional  losses. 

XX.  NUMERICAL  METHOD 

Because  of  the  length  (about  75  pages)  the  code  listing  is  not 
given.  Appendix  A  lists  the  comments  from  the  beginning  of  the  code. 
This  is  a  complete  description  of  the  input  options  and  the  notation 
actually  used  in  the  code  for  the  input  variables.  A  copy  of  the  code 
is  available  from  the  author  on  request. 

The  ordinary  differential  equations  derived  above  are  solved  using 
EPISODE. ^  This  is  a  robust  and  efficient  code  for  the  solution  of 
ordinary  differential  equations.  The  procedure  is  discussed  in  the 
previous  reports. 


XXI.  GE  30MM  GUN  FIXTURE 

As  an  example  of  the  procedure  for  using  the  gun  code,  a  fixture 
developed  by  GE  will  be  modeled.  The  purpose  of  this  specific  gun 
design  was  to  show  that  a  liquid  propellant  regenerative  gun  can  operate 
in  the  high  performance  regime  required  for  tank  cannons.  The  fixture 
has  so  far  only  been  fired  with  a  1/2  charge,  since  problems  arose 
during  the  test  sequence.^  A  schematic  of  the  fixture  is  shown  in 
Figure  2 .  Figure  3  shows  the  experimental  chamber  and  liquid  pressures , 
and  Figure  4  the  piston  travel  (Round  4) .  These 
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Figure  3.  Experimental  Chamber  Pressure  (line)  and  Liquid  Pressure 
(dot)  -  Round  4. 


time  (ms) 


Figure  4.  Experimental  Piston  Travel  -  Round  4 


profiles  will  be  analyzed  in  detail  below.  All  the  data  discussed  in 
this  report  has  been  filtered  to  remove  high  frequency  oscillations. 

12 

In  a  previous  report  test  data  from  an  in-line  regenerative  gun 

was  analyzed  using  an  inverse  code.  Values  of  the  discharge  coefficient 

were  derived.  The  discharge  coefficient  started  small  and  took  several 

milliseconds  to  rise  to  a  value  near  one.  A  similar  analysis  was 

performed  for  the  reverse  gun.  The  vent  area  after  shot  start  was 

computed  based  on  the  experimental  projectile  and  piston  travels.  The 

vent  area  before  shot  start  (through  the  tail  seal)  was  chosen  on  the 

basis  of  agreement  with  the  model.  The  liquid  pressure  measurement 

breaks  down  before  the  end  of  the  firing  cycle,  so  the  liquid  pressure 

is  approximated  as  the  chamber  pressure  times  the  hydraulic  difference. 

In  this  case,  where  the  vent  area  varies  rapidly,  it  is  more  informative 

to  plot  the  discharge  coefficient  times  the  vent  area  (effective 
18 

area) .  Figure  5  compares  the  calculated  actual  area  with  the  derived 
effective  area.  After  shot  start,  the  discharge  coefficient  increases 


Figure  5.  Actual  Vent  Area  (line).  Effective  Area  (dot).  Round  4. 
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much  more  rapidly  chan  for  the  in-line  gun.  But  after  the  vents  are 
fully  open,  the  discharge  coefficient  (effective  area)  decreases  before 
increasing  close  to  the  actual  area.  However,  this  decrease  occurs  at 
the  point  where  the  experimental  chamber  pressure  is  increasing  very 
rapidly  (see  below).  As  a  result,  the  pressures  will  not  have  time  to 
equilibrate.  The  pressure  measured  near  the  wall  of  the  chamber  will 
not  be  equal  to  the  actual  pressure  at  the  outlet  of  the  injectors,  and 

the  liquid  pressure  will  not  equal  the  gun  tube  pressure  times  the 

hydraulic  difference.  Because  of  these  effects,  initially  an  attempt  is 
made  to  model  the  data  assuming  a  large  constant  discharge  coefficient. 

A.  DESCRIPTION  OF  THE  INPUT 

Below  is  a  description  of  the  input  for  the  gun  code.  This 
follows  the  format  in  Appendix  A. 

The  first  line  is  merely  a  label  for  the  problem.  It  lists  the 

filename  of  the  input  job  stream  and  a  brief  description  of  the  problem. 

The  initial  offset  of  the  projectile  and  the  total  distance 
traveled  by  the  projectile  before  muzzle  exit  are  given.  The  offset  is 
chosen  as  the  distance  from  the  start  of  the  short  barrel  to  the  start 
of  the  projectile.  The  diameter  of  the  gun  tube  is  given.  The 
projectile  and  piston  weights  are  entered  (measured). 

The  initial  volumes  of  the  reservoir  and  chamber  are  given.  The 
reservoir  volume  is  computed  from  the  drawings.  The  volume  in  the 
injection  ports  is  ignored.  The  initial  chamber  volume  is  more 
complicated  because  of  the  slots  in  the  stub  case  holder  and  the  short 
barrel.  A  measurement  was  made  of  the  volume  from  the  stub  case  holder 
back  by  filling  the  assembly  with  water. ^  The  rest  of  the  volume  is 
computed.  The  initial  gun  tube  volume  (computed  as  the  projectile 
offset  times  the  tube  area)  is  subtracted  from  the  total  gas  volume. 
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The  areas  of  the  piston  on  the  reservoir  and  chamber  sides  are 
given,  as  well  as  the  grease  dyke  area.  The  gun  tube  area  is  computed 
from  the  gun  tube  diameter. 

The  vent  area  in  the  piston  is  zero.  The  area  between  the  chamber 
and  the  gun  tube  is  estimated  from  the  engineering  drawings. 

The  areas  of  the  three  sets  of  holes  are  given.  Each  set  consists 
of  twelve  holes.  The  length  of  the  holes  is  only  necessary  if  the 
transient  injection  model  is  used. 

The  offset  is  given  from  the  initial  postion  of  the  piston  to  the 
center  of  the  piston  side  injection  holes.  This  is  just  the  maximum 
piston  travel  minus  the  radius  of  the  holes.  As  the  piston  covers  these 
holes,  the  injection  area  will  be  reduced. 

The  offset  from  the  start  of  the  tube  to  the  center  of  both  sets 
of  tube  side  holes  are  given.  The  distance  to  the  first  set  of  holes 
is  chosen  so  that  the  vents  are  just  slightly  open.  The  value  was 
chosen  to  match  the  early  piston  travel.  The  distance  was  the  first  set 
of  holes  to  the  second  set  was  read  off  the  engineering  drawings. 

The  piston  resistance  is  taken  to  be  uniformly  zero.  The 
discharge  coefficient  is  set  at  a  constant  0.95.  This  agrees  well  with 
higher  dimensional  transient  models.^ 

The  flow  into  the  gun  tube  is  normally  modeled  as  steady  state 
Bernoulli  flow  (FLUX1) .  The  flow  into  the  gun  tube  is  steady  state 
isentropic  flow  (FLUX1) . 

The  nominal  shot  start  pressure  is  70  MPa,  based  on  the  strength 
of  the  break  link.  The  experimental  data  shows  that  the  projectile 
moves  slightly  sooner,  and  a  value  of  65  MPa  is  used.  After  shot  start 
the  resistance  pressure  is  taken  as  5  MPa.  This  is  a  typical  value  from 
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solid  propellant  gun  modeling. 

.  Next  the  physical  properties  of  the  propellant  HAN1845  are  given. 

The  density  at  atmospheric  pressure  has  been  measured  and  the  bulk 

modulus  has  been  fitted. ^  The  chemical  energy,  the  ratios  of  specific 

heats,  the  molecular  weight  of  the  propellant  final  produces  gas,  and 

o 

the  covolume  have  been  computed  using  BLAKE.  The  surface  tension  and 

O 

the  kinematic  viscosity  of  the  liquid  are  not  actually  required,  but 
are  used  to  compute  the  Weber  number,  the  Reynold's  number,  and  the 
third  Lagrange  number. 

The  liquid  is  pre-presurized  to  about  1  MPa.  The  gas  is  initially 
at  one  atmosphere. 

As  a  first  approximation,  the  liquid  is  assumed  to  combust 
instantaneously  as  it  enters  the  combustion  chamber  (DR0P1) .  Later  the 
effects  of  accumulation  will  be  studied. 

The  primer  is  assumed  to  be  injected  in  the  form  of  hot  gas 
(PRIM3) .  There  are  3.0  grams  of  IMR  powder  in  the  booster.  Experiments 
with  just  the  primer  (water  and  dry  shots)  show  an  injection  time  of  1.5 
ms.  The  heat  loss  was  chosen  to  obtain  a  pressure  due  to  the  primer 
between  the  water  and  dry  shot  results. 

Heat  loss  to  the  gun  tube  walls  is  calculated.  The  gun  tube  wall 
temperature  is  taken  to  be  300K.  The  heat  loss  algorithm  in  the  code  is 
not  adjusted. 

The  air  shock  in  front  of  the  projectile  is  calculated.  The  air 
in  the  tube  is  initially  assumed  to  be  at  1  atm.  and  300K.  The  values  ' 
of  the  molecular  weight  and  ratio  of  specific  heats  are  taken  at  these 
conditions. 

The  standard  Lagrange  pressure  distribution  in  the  gun  tube  is 
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assumed. 

The  code  will  print  out  results  every  0.05  ms.  (TINC) .  Because 
the  code  must  often  change  the  time  step,  it  is  more  efficient  to 
restrict  the  maximum  time  step  (HTOF) .  The  error  controls  EPS  and  SREC 
are  given  typical  values. 

The  integration  method  flag  HF  is  set  to  22  (backwards 
differentiation  formulas  with  an  internally  computed  Jacobian) .  KWR1TE 
is  set  to  zero  to  eliminate  diagnostic  messages.  A  time  limit  TMAX  is 
set.  If  the  code  takes  longer  than  TMAX  seconds  to  execute,  the  code 
will  stop  and  write  the  usual  summary  pages  so  all  the  information  will 
not  be  lost.  The  code  is  only  to  be  integrated  once  (REP1)  and  the 
chamber  pressure  will  be  computed  normally  (CHAMl) 

An  earlier  version  of  these  computations  was  presented  at  a  JANNAF 

OQ  * 

meeting.  The  input  used  at  that  time  was  somewhat  different.  The 
hardware  for  this  fixture  recently  arrived  at  BEIL,  and  was  somewhat 
different  than  the  engineering  drawings.  Also,  there  was  an  error  in 
the  earlier  calculation  of  the  initial  gas  volume.  The  primary 
difference  between  this  and  the  previous  calculations  is  the  smaller 
initial  gas  volume.  Also,  the  grease  dyke  is  smaller  than  in  the 
previous  calculations,  which  causes  a  larger  hydraulic  difference. 

B .  PRELIMINARY  RUN 

A  preliminary  run  is  made  using  the  experimental  chamber  pressure 

(CHAM3) .  For  numerical  reasons,  a  spline  fit  of  the  pressure  is  used 

instead  of  the  raw  data  (see  Figure  6).  The  initial  projectile  offset 

is  adjusted  to  match  the  early  piston  travel.  Figure  7  shows  the 

resulting  piston  travel.  The  agreement  is  quite  good.  Unlike  the  in- 
o 

line  gun,  a  constant  large  discharge  coefficient  is  a  good 
approximation. 
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Figure  6.  Experimental  Chamber  Pressure  -  Round  4  (line) 
Spline  Fit  to  the  Chamber  Pressure  (dot). 


Figure  7.  Experimental  Piston  Travel  -  Round  4  (line). 
Model  with  Experimental  Chamber  Pressure  (dot). 


The  effect  of  recoil  is  a  possible  problem.  The  piston  travel  was 
measured  in  a  relative  sense.  The  data  is  normalized  to  the  maximum 
piston  travel.  However,  the  gun  recoil  is  not  measured.  To  estimate 
the  size  of  this  effect,  the  recoil  is  computed  based  on  the  forces 
acting  on  the  gun.  The  weight  of  the  gun  is  estimated  from  its 
volume.  Since  the  recoil  mechanism  has  not  been  studied,  free  recoil 
is  assumed.  The  computed  recoil  will  then  be  larger  than  the  actual 
recoil.  Figure  8  shows  the  computed  recoil,  compared  to  the 
experimental  and  model  piston  travel.  Recoil  does  not  begin  to  be 
important  until  near  the  end  of  the  piston  stroke.  The  model  shows  a 
very  slow  approach  by  the  piston  to  the  end  of  the  reservoir.  The 
experiment  instead  shows  a  small  reversal.  The  size  of  this  reversal  is 
in  reasonable  agreement  with  the  computed  recoil.  The  later  and  much 
faster  reversal  cannot  be  due  to  recoil,  and  must  instead  be  actual 

piston  motion.  The  final  result  is  that  the  experimental  piston  travel 

is  slightly  too  large.  The  effect  is  small,  and  no  attempt  is  made  to 

correct  for  the  effect  of  recoil.  Since  the  mount  restricts  the  recoil, 

the  actual  error  due  to  recoil  is  very  small. 


Figure  8.  Experimental  Piston  Travel  -  Round  4  (line). 
Computed  Gun  Recoil  (dot). 
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The  muzzle  velocity  predicted  by  the  code  is  more  than  10%  high. 
The  Lagrange  pressure  distribution  is  derived  under  the  assumption  that 
the  pressure  equilibrates  instantly  in  the  gun  tube.  In  the  case 
considered  here,  the  pressure  increase  is  too  rapid,  and  the  pressure  at 
the  base  of  the  projectile  will  be  lower  than  that  predicted  by  the 
Lagrange  distribution  (see  below). 

In  earlier  work  on  this  fixture  by  Bulman  and  Maher  ,  a  low 
effective  area  (discharge  coefficient)  was  blamed  for  the  delayed 
propellant  ignition.  Poor  flow  through  the  "Y"  shaped  passages  was 
considered  to  be  the  major  problem,  rather  than  liquid  accumulation. 

This  is  opposite  the  conclusion  reached  in  this  paper.  As  a  partial 
check,  the  code  was  run  still  using  the  experimental  chamber  pressure 
but  assuming  a  discharge  coefficient  of  0.85  and  0.75  (Figure  9).  The 
initial  projectile  offset  was  adjusted  slightly  for  the  two  new  cases  to 
match  the  early  piston  travel.  There  is  not  a  great  deal  of  difference 
in  the  profiles.  As  mentioned  above,  the  experimental  piston  travel  is 


Figure  9.  Experimental  Piston  Travel  -  Round  4  (line). 
Model  with  Experimental  Chamber  Pressure.  -  0.95  (dot). 
Cq  -  0.85  (dash).  -  0.75  (dot -dash). 
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not  known  completely  accurately.  Also,  the  pressure  rise  is  too  rapid 
for  the  system  to  equilibrate.  The  pressure  at  the  injection  holes 
outlet  will  probably  be  higher  in  reality  than  in  the  model,  leading  to 
a  lower  apparent  discharge  coefficient.  So  the  value  of  0.95  is  still 
thought  to  be  the  best  approximation  to  the  system.  However,  there  may 
be  minor  flow  problems  in  the  "Y"  shaped  passages. 

C.  MODELING  -  ROUNDS  4-10 

The  model  is  used  to  compute  the  pressure,  assuming  instantaneous 
burning.  Figures  10  and  11  show  the  results.  To  make  comparisons  more 
easily,  the  model  curves  have  been  shifted  uniformly  in  time  so  the 
chamber  pressures  match  at  75  MPa.  The  agreement  is  very  poor.  This 
run  shows  that  liquid  accumulation  plays  a  major  role  in  the  fixture. 
However,  this  model  does  reproduce  the  first  part  of  the  rapid  chamber 
pressure  rise. 

The  droplet  option  was  introduced  next.  The  drop  size  versus 
projectile  travel  table  was  adjusted  by  trial  and  error  until  reasonable 
agreement  with  the  chamber  pressure  was  obtained.  The  change  from  the 
DR0P1  option  to  DR0P4  is  the  only  change  made  in  the  code.  The  results 
(Figures  12-13)  no  longer  require  centering,  since  the  delay  in  the 
experimental  pressure  rise  is  now  being  modeled.  The  agreement  with 
the  experimental  data  is  quite  good.  Figure  14  shows  the  droplet 
diameter,  and  Figure  15  the  calculated  liquid  accumulation  in  the 
chamber/  gun  tube. 

Appendix  B  shows  the  input  for  this  run.  The  numbers  and  labels 
at  the  left  are  read  by  the  code.  The  comments  at  the  right  are  for 
user  identification.  The  corresponding  output  is  also  given.  The  code 
was  run  on  the  BRL  Cray-2.  The  input  is  echoed  back,  along  with  other 
calculated  initial  values.  Once  the  integration  begins,  output  is  given 
every  0.05  ms.  The  pressures,  piston  travel  and  velocity,  and 
projectile  travel  and  velocity  are  given.  At  the  end  of  the 
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13.  Experimental  Piston  Travel  -  Round 
Droplet  burning  (dot). 
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Figure  15.  Experimental  Chamber  Pressure  -  Round  4  (line). 

Mass  in  the  Chamber/  Gun  Tube  (dot).  Liquid  Accumulation  (dash). 
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integration,  a  summary  page  is  given.  A  second  page  shows  when  various 
maximum  values  are  achieved.  There  are  many  other  variables  of  interest 
in  the  simulation.  As  the  integration  proceeds,  these  are  written  onto 
other  files,  and  can  be  appended  to  the  output.  In  Appendix  B,  only  the 
minimal  output  is  given. 

Now  consider  the  profiles  in  more  detail.  The  primer  causes  the 
initial  pressure  rise.  Liquid  is  then  injected  slowly  around  the  tail 
seal  on  the  projectile.  According  to  the  model,  only  about  1/5  of  the 
liquid  injected  before  shot  start  actually  combusts.  At  about  65  MPa, 
the  shot  start  link  releases.  The  pressure  now  decreases.  The  rapid 
projectile  motion  opens  the  vents  completely.  It  also  causes  the  hot 
gas  to  cool  by  isentropic  expansion.  The  cold  liquid  propellant 
suddenly  is  injected  into  the  cooler  gas  and  quenches  the  combustion. 
Combustion  is  almost  negligibly  small.  Since  liquid  is  being  injected 
rapidly,  the  liquid  accumulation  becomes  large  (about  1/3  of  the 
original  charge).  Eventually  the  liquid  reignites,  causing  a  very  rapid 
pressure  rise.  The  measured  chamber  pressure  increases  90  MPa  in  0.14 
ms.  The  pressure  increase  in  the  gun  tube  itself  is  not  measured,  and 
may  be  even  more  steep.  Combustion  is  so  rapid  compared  to  the  speed  of  . 
sound  that  large  pressure  gradients  will  be  created.  The  model 
assumptions  of  uniform  pressure  in  the  chamber  and  a  Lagrange  pressure 
distribution  in  the  gun  tube  are  not  accurate. 

Only  a  small  amount  of  the  accumulation  bums  off  at  this  time. 

The  experiment  does  show  a  higher  pressure  rise  than  the  model,  so  more 
liquid  will  be  burned  off  in  practice.  However,  a  large  part  of  the 
pressure  rise  is  injection  driven.  The  local  pressure  rises  so  rapidly 
that  the  liquid  pressure  does  not  have  time  to  equilibrate.  The  local 
gun  tube  pressure  than  can  become  as  large  as  the  liquid  pressure  (see 
Figure  3)  and  injection  is  shut  off.  The  chamber  pressure  than  falls 
steeply,  since  the  combustion  is  largely  injection  driven.  The  pressure 
shows  low  frequency  oscilations  for  the  rest  of  stroke.  Since  the  model 
assumes  that  pressures  equilibrate  instantly,  it  does  not  reproduce  this 
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behavior . 


The  plscon  Chen  moves  over  Che  injecCion  holes.  This  causes  a 
rapid  deceleracion  of  Che  piscon.  Only  near  che  end  of  Che  injeccion 
process  does  mosC  of  che  accumulaCion  bum  off.  EvenCually  Che  hoc  gas 
propagaCes  back  inCo  che  reservoir.  In  Chis  case,  Chere  is  very  liCCle 
liquid  lefc,  and  Che  reversal  is  very  mild,  as  shown  in  Figure  8  aC 
abouC  7  ms. 

Two  ocher  shoes  (Rounds  7  and  8)  gave  resulcs  similar  Co  Round  4. 
For  Round  6,  Che  break  link  released  premaCurely.  Due  Co  Che  lower 
pressure,  che  liquid  never  reigniCed  (flameouc). 

Round  10  is  similar  Co  Round  4  (see  Figure  16).  However,  chere 
was  an  even  longer  delay  before  che  liquid  reigniCed.  The  resulcing 
chamber  pressure  rise  was  very  rapid  (150  MPa  in  0.04  ms).  There  was  a 
piscon  reversal  during  Chis  sceep  pressure  rise,  alchough  blow-ouC  plugs 
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in  the  reservoir  prevented  damage  to  the  system.  Since  most  of  the 
combustion  takes  place  in  the  tube,  the  gun  tube  pressure  became  much 
higher  than  the  chamber  pressure  at  the  piston,  and  even  higher  than  the 
liquid  pressure.  Hot  gas  then  easily  enters  the  reservoir  and  causes 
the  piston  reversal. 

The  major  problem  with  this  fixture  is  the  large  liquid 
accumulation.  This  can  lead  to  flameout,  or  to  a  very  rapid  pressure 
rise  if  the  liquid  reignites.  The  pressure  rise  will  be  rapid  even 
without  accumulation  (see  Figure  10) ,  because  the  initial  gas  volume  is 
small.  However,  a  rapid  local  pressure  rise  will  shut  off  the 
injection.  The  presence  of  accumulated  liquid  allows  the  pressure  to 
keep  increasing  after  injection  has  stopped,  causing  the  back  flow  and 
the  piston  reversal. 

D.  MODELING  -  ROUNDS  25-28 


An  earlier  fixture,  also  without  injection  through  the  piston,  was 
more  successful.  The  present  fixture  was  modified  to  bring  it  closer 
to  the  previous  desitn.  The  injection  area  was  cut  in  half.  The 
forward  and  backward  holes  were  staggered,  instead  of  meeting  at  a  "Y" . 
The  projectile  mass  was  increased  to  350  grams.  A  different  break  link 
was  used,  with  a  shot  start  of  45  MPa.  The  break  link  also  stretched, 
leading  to  a  proportionally  larger  initial  injection  area.  Slightly 
less  primer  was  used. 

Figures  17  and  18  show  the  experimental  data  for  Round  25.  There 
is  no  longer  a  delay  in  the  pressure  rise  after  shot  start.  It  is 
harder  to  compare  the  piston  travel  because  recoil  is  more  important  for 
this  design.  Because  of  the  heavier  projectile  and  hence  lower 
velocity,  higher  pressures  are  exerted  for  a  longer  period  of  time. 

As  before,  a  preliminary  run  was  made  using  the  experimental 
chamber  pressure  (CHAM2  option) .  The  agreement  with  the  experimental 
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muzzle  velocity  is  now  very  good.  The  Lagrange  pressure  distribution  is 
adequate  if  the  pressure  rise  is  not  too  rapid. 

Using  the  droplet  option,  this  case  can  be  modeled  more  simply 
than  Round  4.  Appendix  C  gives  the  input  and  the  summary  page  for  the 
output.  Combustion  is  still  inefficient  at  first.  But  when  the 
projectile  starts  to  move,  the  liquid  is  injected  slowly  enough  so  that 
relatively  little  accumulation  occurs  (Figures  19-22).  The  agreement  in 
piston  travel  is  good.  If  the  experimental  piston  travel  could  be 
corrected  for  recoil,  the  agreement  would  probably  improve.  This  cannot 
presently  be  done  because  the  effect  of  the  mount  on  the  actual  recoil 
is  unknown . 

The  experiment  still  shows  a  rapid  pressure  rise  just  before  3.0 
ms.  This  occurs  right  after  the  projectile  clears  the  injection  holes. 
Since  the  projectile  is  much  heavier  and  the  shot  start  is  lower,  there 
is  more  of  a  delay  between  shot  start  and  the  opening  of  the  injectors. 
The  pressure  Increases  rapidly  for  a  short  time  and  then  levels  off. 

The  pressure  rise  is  as  before  too  rapid  for  the  liquid  pressure  to 
equilibrate.  The  rapid  local  pressure  rise  slows  down  dramatically  the 
injection  rate.  The  accumulation  is  not  burning  off  at  this  point,  so 
when  the  injection  rate  becomes  slow  enough,  the  pressure  stops 
increasing.  This  occurs  at  least  once  more  before  peak  pressure.  The 
model,  since  it  assumes  that  the  pressures  equilibrate,  does  not 
demonstrate  this  behavior. 

Rounds  26  and  27  were  very  similar  to  Round  25.  However,  Round  28 
showed  a  piston  reversal.  The  reversal  occured  soon  after  the  peak 
pressure,  as  the  piston  was  just  starting  to  block  off  the  injection 
holes.  The  piston  actually  starts  to  decelerate  compared  with  the  model 
right  after  the  projectile  clears  the  injection  ports.  At  this  point 
the  chamber  pressure  goes  up  rapidly  compared  to  Round  25  (40  MPa  in  0.1 
ms).  The  recorded  chamber  pressure  is  essentially  equal  to  the  recorded 
liquid  pressure,  and  the  gun  tube  pressure  should  be  higher  than  the 
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Piston  travel  I  cm) 


Figure  20.  Experimental  Piston  Travel  -  Round  25  (line). 
Model  -  Droplet  Burning  (dot) . 
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Figure  21.  Model  -  Droplet  Diameter  -  Round  25. 


Figure  22.  Experimental  Chamber  Pressure  -  Round  25  (line). 

Mass  in  the  Chamber/  Gun  tube  (dot).  Liquid  Accumulation  (dash) 


chamber  pressure  (see  Figure  23).  Hot  gas  then  could  start  to  back  up 
into  the  injection  ports.  The  pressure  rise  is  much  less  dramatic  than 
for  Round  10,  and  the  corresponding  piston  reversal  is  more  gradual. 


The  initial  pressure  rise  (before  shot  start)  is  more  rapid  here 
than  for  Round  23,  implying  more  efficient  combustion.  The  same 
injection  driven  pressure  rise  occurs  when  the  projectile  opens  the 
injection  ports.  However,  in  this  case  a  small  amount  of  the 
accumulated  liquid  also  burns ,  pushing  the  pressure  high  enough  to  cause 
the  gas  backup. 

There  is  another  important  point  here.  Since  there  is  no  longer  a 
quenching  period  before  the  rapid  pressure  rise,  the  volume  is  only 
slightly  larger  than  the  initial  very  small  gas  volume  of  50  cm^  when 
the  pressure  rise  occurs.  It  then  takes  less  combustion  to  cause  a 
given  pressure  rise.  It  is  easier  to  cause  a  rapid  pressure  rise  and 
gas  backup. 


Figure  23.  Experimental  Chamber  Pressure  (line)  and 
Liquid  Pressure  (dot)  -  Round  28. 


The  previous  gun  fixture  did  not  exhibit  the  sane  problems  with 

I  O 

piston  reversal.  However,  the  earlier  fixture  had  a  larger  hydraulic 
difference.  The  liquid  pressure  would  be  higher  compared  to  the  gun 
tube  pressure,  making  a  reversal  less  likely. 

E.  DISCUSSION 

The  experimental  data  can  be  reproduced  accurately  using  the 
model.  Assuming  a  large  constant  discharge  coefficient  is  a  good 
approximation,  unlike  the  in-line  case.  However,  liquid  accumulation  is 
very  important,  especially  for  the  larger  injection  area. 

The  major  problem  with  this  design  is  the  rapid  injection  of 
liquid  into  a  small  volume  as  the  projectile  first  starts  to  move.  This 
may  totally  quench  the  combustion,  causing  a  flameout,  or  it  may 
temporarily  quench  the  combustion.  When  the  liquid  reignites,  the  very 
rapid  pressure  increase  leads  to  non-uniform  pressure  and  a  possible 
dramatic  piston  reversal.  Even  if  the  combustion  does  not  quench,  the 
rapid  injection  leads  to  a  rapid  pressure  increase.  Again,  the  non- 
uniform  pressure  allows  hot  gas  to  penetrate  into  the  reservoir,  in  this 
case  causing  a  much  milder  piston  reversal. 

The  fixture  is  designed  to  reach  high  muzzle  velocities.  This 
implies  that  liquid  must  be  injected  rapidly.  As  the  fixture  moves 
closer  to  the  high  performance  regime,  the  problems  can  be  expected  to 
get  worse. 

One  previous  suggestion  is  a  multi-angle  injector.^  The  liquid 
would  be  introduced  somewhat  more  gradually  and  would  be  spread  out  more 
along  the  gun  tube.  The  fixture  would  be  more  complicated.  And  since 
the  projectile  moves  rapidly  at  shot  start,  all  the  injection  ports 
would  be  opened  very  quickly. 

A  simpler  idea  is  to  change  the  design  of  the  booster  and/or  the 
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tail  seal.  At  present  the  hot  igniter  gases  are  aimed  directly  at  the 
piston.  Meanwhile,  cold  liquid  is  injected  around  the  tail  seal  and 
directed  along  the  wall  of  the  short  barrel.  The  model  shows  that  very 
little  of  this  liquid  burns.  Since  it  will  be  cooled  by  contact  with 
the  wall,  this  is  not  surprising.  Some  of  the  holes  in  the  booster 
could  be  redirected  toward  the  wall  of  the  short  barrel.  This  would 
tend  to  ignite  the  liquid  and  prevent  the  early  accumulation.  In 
addition,  the  tail  seal  could  be  redesigned  to  direct  the  liquid  toward 
the  center  of  the  gun  tube  rather  than  along  the  walls.  If  the  liquid 
can  be  ignited  early,  the  combustion  is  not  likely  to  be  quenched  when 
the  main  injection  process  begins. 

Even  if  the  quenching  phase  is  eliminated,  this  may  not  completely 
solve  the  problem,  as  Round  28  shows.  There  can  be  a  rapid  pressure 
rise  even  with  very  little  accumulation.  The  initial  gas  volume  could 
be  made  larger,  which  would  make  it  harder  for  the  pressure  to  rise 
rapidly.  Also,  the  piston  could  be  redesigned  to  obtain  a  larger 
hydraulic  difference.  If  the  liquid  pressure  is  designed  to  be  higher 
relative  to  the  chamber  pressure,  gas  backup  into  the  reservoir  becomes 
less  likely. 

As  a  last  resort,  some  of  the  liquid  could  be  injected  through  the 
piston.  This  was  done  for  most  of  the  earlier  tests  conducted  with  this 
type  of  fixture.  With  propellant  burning  directly  outside  the 
reservoir,  the  liquid  pressure  will  stay  higher.  Also,  if  the  fixture 
is  designed  so  there  is  gas  flow  out  of  the  chamber  into  the  tube,  this 
will  tend  to  break  up  liquid  concentrations  near  the  gun  tube  injection 
holes.  However,  this  would  complicate  the  fixture  considerably. 

XXII.  CONCLUSIONS 

A  description  is  given  of  a  lumped  parameter  model  for  reverse 
annular  regenerative  liquid  propellant  guns.  A  complete  discussion  of 
the  governing  equations  and  the  different  options  is  given.  The  code 


56 


has  been  used  Co  test  Che  effecc  of  various  assumptions,  so  addicional 
opcions  will  be  added  as  needed. 

The  model  resulcs  are  compared  vich  experimenCal  daCa  from  a  30mm 
gun  fixCure.  A  discussion  is  given  of  Che  procedure  for  choosing  a  sec 
of  input  parameCers.  Assuming  a  large  sCeady  discharge  coefficienC  is  a 
good  approximaCion,  unlike  Che  case  of  in-line  guns.  By  assuming 
instanCaneous  burning,  Che  overall  performance  of  Che  gun  can  be 
predicced  wich  reasonable  accuracy. 

Liquid  accumulaCion  is  very  imporcanc  for  Chis  fixcure.  This 
cannoC  be  predicced,  since  Che  flow  fields  and  combusCion  are  noC  well 
underscood,  buc  Che  code  can  be  adjusted  Co  macch  a  given  experimenc. 

The  code  can  Chen  be  used  Co  sCudy  decails  noC  available  from  Che 
experimenc . 

This  parcicular  fixcure  exhibits  difficulties.  Just  after  shot 
start  there  can  be  excessive  liquid  accumulaCion,  quenching  the  flame. 
After  Che  maximum  pressure  there  can  be  piston  reversal.  Neither  effecc 
can  be  predicted  by  the  lumped  parameter  code,  but  the  model  results 
help  explore  different  hypotheses  for  this  behavior. 
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GLOSSARY 


o 

Area  of  the  piston  -  liquid  side,  cm  . 

2 

A2  Area  of  the  piston  •  chamber  side,  cm  . 

2 

A^  Area  of  the  gun  tube,  cm  . 

2 

Ag  Area  of  the  grease  dyke  around  the  piston,  cm  . 

2 

A^2  Area  of  the  vent  into  the  chamber  from  the  reservoir,  cm  . 

2 

Aj^  Area  of  the  vent  into  the  gun  tube  from  the  reservoir,  cm 
b  Covolume ,  cm  /g . 

c^  Speed  of  sound  in  the  liquid  in  the  reservoir,  cm/s. 

C2  Speed  of  sound  in  the  mixture  in  the  chamber,  cm/s. 

c^  Speed  of  sound  in  the  mixture  in  the  gun  tube,  cm/s. 

c^2  Speed  of  sound  in  the  liquid  in  the  chamber,  cm/s. 

CL4  Speed  of  sound  in  the  liquid  in  the  gun  tube,  cm/s. 

Cq2  Speed  of  sound  in  the  gas  in  the  chamber,  cm/s. 

CG4  Average  speed  of  sound  in  the  gas  in  the  gun  tube,  cm/s. 

Cp  Specific  heat  at  constant  pressure,  j/g-K. 

cv  Specific  heat  at  constant  volume,  j/g-K. 
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c 

d, 


el 

e2 


e4 

EL1 

EL2 

EL4 

EG2 

EG4 

EP 

«1 A 
^4 

“pj 

EH4 


Discharge  coefficient  for  the  mass  flux  from  the  reservoir. 
Diameter  of  the  gun  tube,  cm. 

Chemical  energy  of  the  liquid,  j/g. 

Internal  energy  of  the  gas  in  the  chamber,  j/g. 

Internal  energy  of  the  gas  in  the  gun  tube,  j/g. 

Total  internal  energy  in  the  liquid  reservoir,  j. 

Total  internal  energy  in  the  liquid  in  the  chamber ,  j . 

Total  internal  energy  in  the  liquid  in  the  gun  tube,  j. 

Total  internal  energy  in  the  gas  in  the  chamber ,  j . 

Total  internal  energy  in  the  gas  in  the  gun  tube,  j. 

Total  internal  energy  in  the  unbumed  primer ,  j  . 

Kinetic  energy  of  the  liquid  in  the  gun  tube,  j. 

Kinetic  energy  of  the  gas  in  the  gun  tube,  j. 

Kinetic  energy  of  the  projectile,  j. 

Kinetic  energy  of  the  piston,  j . 

Total  heat  loss  to  the  gun  tube  walls ,  j . 

Total  energy  loss  due  to  projectile  friction,  j. 
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EF. 


ps 


So 

hLl 

hL2 

hL4 

hG2 

hG4 

*S, 

K 


m12 

“14 

m24 

V 

m2 


Total  energy  loss  due  to  piston  friction,  j . 

Total  energy  of  the  system,  j. 

7  2 

Conversion  constant  -  10  g/MPa-cm-s  . 

Liquid  enthalpy  in  the  reservoir,  j/g. 

Liquid  enthalpy  in  the  chamber,  j/g. 

Liquid  enthalpy  in  the  gun  tube,  j/g. 

Gas  enthalpy  in  the  chamber,  j/g. 

Gas  enthalpy  in  the  gun  tube,  j/g. 

Heat  transfer  coefficient  to  the  gun  tube  walls,  j/cm  -K-s. 
Bulk  modulus,  MPa. 

Bulk  modulus  at  zero  pressure,  MPa. 

Derivative  of  the  bulk  modulus . 

Mass  flux  from  the  reservoir  into  the  chamber,  g/s. 

Mass  flux  from  the  reservoir  into  the  gun  tube,  g/s. 

Mass  flux  from  the  chamber  into  the  gun  tube,  g/s. 

Mass  flux  of  the  primer  into  the  gun  tube,  g/s. 

Chamber  rate  of  production  of  gas  by  combustion,  g/s. 
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“4 

mL2 

“L4 

®G4 

m2 

M4 

ML2 

ML4 

“G2 

«G4 

"p 

Mj 

*pj 

Hps 


Gun  tube  rate  of  production  of  gas  by  combustion,  g/s. 
Total  rate  of  production  of  liquid  in  the  chamber,  g/s. 

Total  rate  of  production  of  liquid  in  the  gun  tube,  g/s. 

Total  rate  of  production  of  gas  in  the  chamber,  g/s. 

Total  rate  of  production  of  gas  in  the  gun  tube,  g/s. 

Total  mass  in  the  reservoir,  g. 

Total  mass  in  the  chamber,  g. 

Total  mass  in  the  gun  tube,  g. 

Liquid  mass  in  the  chamber,  g. 

Liquid  mass  in  the  gun  tube,  g. 

Gas  mass  in  the  chamber,  g. 

Gas  mass  in  the  gun  tube,  g. 

Unburned  primer  mass,  g. 

Total  propellant  mass  in  the  system,  g. 

Mass  of  the  projectile,  g. 

Mass  of  the  piston,  g. 

Molecular  weight  of  the  gas,  g/mole. 
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Pi 

P2 

P4 

PL 

PR 

PPj 

Pps 


R, 


*u 


*PJ 


ps 


Pressure  in  the  liquid  reservoir,  MPa. 

Pressure  in  the  chamber,  MPa. 

Space  mean  pressure  in  the  gun  tube,  MPa. 

Pressure  at  the  breech  end  of  the  gun  tube,  MPa. 

Pressure  at  the  base  of  the  projectile,  MPa. 

Projectile  resistance  pressure,  MPa. 

Piston  resistance  pressure,  MPa. 

Air  shock  pressure  ahead  of  the  projectile,  MPa. 

Heat  loss  to  the  gun  tube  walls,  j/cm  -s. 

Specific  gas  constant,  j/g-K. 

Universal  gas  constant  -  8.318  j/mole-K. 

Projectile  travel,  cm. 

Piston  travel,  cm. 

Chamber  mass  production  rate  times  enthalpy  change  ,  j/s. 
Gun  tube  mass  production  rate  times  enthalpy  change  ,  j/s. 
Time ,  s . 

Temperature  in  the  chamber,  K. 
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Average  temperature  in  the  gun  tube,  K. 

Tw  Temperature  of  the  gun  tube  walls,  K. 

v12  Injection  velocity  of  the  liquid  into  the  chamber,  cm/s. 

Vj^  Injection  velocity  of  the  liquid  into  the  gun  tube,  cm/s. 

V24  Injection  velocity  from  the  chamber  into  the  gun  tube,  cm/s. 
Vpj  Velocity  of  the  projectile,  cm/s. 

Vps  Velocity  of  the  piston,  cm/s. 

3 

V^  Volume  of  the  liquid  reservoir,  cm  . 

3 

V2  Volume  of  the  chamber,  cm  . 

O 

V^  Volume  of  the  gun  tube  behind  the  projectile,  cm  . 

VL2  Volume  of  the  liquid  in  the  chamber,  cm  . 

3 

VL4  Volume  of  the  liquid  in  the  gun  tube,  cm  . 

3 

Vq2  Volume  of  the  gas  in  the  chamber,  cm  . 

3 

VG4  Volume  of  the  gas  in  the  gun  tube,  cm  . 

xR  Distance  from  the  tube  end  to  the  projectile  base,  cm. 

7  Ratio  of  specific  heats. 

«2  Porosity  of  the  mixture  in  the  chamber. 
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e4 


A 

^o 


n 


p2 


p  4 


p12 


pl A 


PG2 

PGU 


Porosity  of  the  mixture  in  the  gun  tube 
Impetus  of  the  propellant,  j/g. 

3 

Liquid  density  at  zero  pressure,  g/cm  . 

3 

Liquid  density  in  the  reservoir,  g/cm  . 
Mixture  density  in  the  chamber,  g/cm  . 

3 

Mixture  density  in  the  gun  tube,  g/cm  . 

3 

Liquid  density  in  the  chamber,  g/cm  . 

3 

Liquid  density  in  the  gun  tube,  g/cm  . 

3 

Gas  density  in  the  chamber,  g/cm  . 

3 

Gas  density  in  the  gun  tube,  g/cm  . 
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APPENDIX  A 
INPUT  OPTIONS 


The  computer  code  is  long  (around  75  pages) ,  so  a  complete 
listing  is  not  given.  Instead,  only  the  description  of  the 
input  options  from  the  beginning  of  the  code  is  given.  A 
complete  listing  of  the  code  (tape,  disk,  or  hardcopy)  may 
be  obtained  from  the  author. 


c  ********************************************* 
c  December  1,  1988. 
c  Description  of  input 
c 

c  file  1:  problem  title 

c  80  alphanumeric  characters 

c  will  be  title  on  graphics  files 
c 

c  file  2:  offset,  travel 

c  offset  =  distance  from  end  of  tube  to  projectile  (cm) 

c  travel  =  total  projectile  travel  (cm) 

c 

c  file  3:d4 

c  d4  =  gun  tube  diameter  (cm) 

c 

c  file  4:  pjwt 

c  pjwt  *  projectile  weight  (g) 

c 

c  file  5:  pswt 

c  pswt  »  piston  weight  (g) 

c 

c  file  6:  vli,v2i 

c  vii  *  initial  liquid  reservoir  volume  (cm**3) 

c  v2i  *  initial  intermediate  chamber  volume  (cm**3) 

c 

c  file  7:  al,a2,agres 

c  al  *  area  of  liquid  side  piston  face  (including  vents)  (cm**2) 

c  a2  *  area  of  chamber  side  piston  face  (including  vents)  (cm**2) 
c  agres  =  area  of  grease  dyke  (cm* *2) . 
c 

c  file  8:  avl2,av24 

c  avl2  =  vent  area  between  reservoir  and  chamber  (cm**2) . 

c  av24  =  vent  area  between  chamber  and  tube  (cm**2) . 

c 

c  file  9:  avl,av4a,av4b 

c  avl  -  area  y  vent  on  liquid  side  (cm**2) 

c  av4a  *  area  first  y  vent  hole  on  tube  side  (cm**2) 

c  av4b  =  area  second  y  vent  hole  on  tube  side  (cm**2) 

c 

c  file  10:  nvo,  pth 

c  nvo  *  no.  of  holes  between  reservoir  and  tube, 
c  total  vent  area  *  avl4  *  nvo. 
c  avl4  -  min (avl,  av4a+av4b) 

c  pth  =  length  of  the  holes  (thickness  of  short  barrel)  (cm) . 
c 
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c  file  H:sol,so4a,so4b 

c  sol  **  piston  travel  to  center  of  vent  avl  (cm)  . 

c  soa4  =  distance  from  start  of  gun  tube  to  center  of  av4a  (cm) 

c  sob4  *  distance  from  start  of  gun  tube  to  center  of  av4b  (cm) 

c 

c  file  12:  tpis  -  piston  friction  resistance  option 
c 

c  if  pisl  then  resistance  is  a  function  of  piston  travel 

c  file  12.1:  npis 

c  npis  »  no.  of  entries  in  piston  resistance  table 

c  file  12.2...:  sfrac(i) ,prs(i) 

c  sfrac(i)  »  fraction  of  piston  travel 

c  will  be  normalized  to  max  piston  travel 

c  prs(i)  -  piston  resistance  (MPa) 

c  piston  resistance  is  found  by  interpolating  table 

c 

c  file  13:  tdis  -  discharge  coefficient  option  switch 
c 

c  if  disl  then  dc  function  of  piston  travel 

c  file  13.1:  ndc 

c  ndc  *  no.  of  table  entries 

c  file  13.2...:  sfrac(i) ,dcf (i) 

c  sfrac(i)  =  fraction  of  piston  travel 

c  will  be  normalized  to  max  piston  travel 

c  dcf(i)  -  discharge  coefficient 

c  dc  is  found  by  interpolating  table 

c  dc  applied  to  flow  from  reservoir  to  chamber  and  to  tube, 

c 

c  file  14:  tflux  -  mass  flux  option  switch  for  orifice 
c 

c  if  fluxl  then  steady  state  formulation 

c  fluid  entrance  velocity  *  0 

c 

c  if  flux2  then  unsteady  mass  flux  formulation 

c  take  time  derivative  of  mass  flux 

c  assume  mass  flux  constant  wrt  x 

c  file  14.1:  fudge 

c  fudge  «  fudge  factor 

c  divide  time  derivative  of  mass  flux  by  fudge 

c 

c  file  15:  tfluxp  -  mass  flux  option  switch  for  tube 
c 

c  if  fluxl  then  steady  state  formulation 

c  generalized  bemoulli  with  energy  difference 

c  isentropic  flow 

c  fluid  entrance  velocity  =  0 

c 

c  file  16:  tproj  -  projectile  resistance  option  switch 
c 

c  if  projl  then  projectile  resistance  function  of  proj  travel 

c  file  16.1:  nproj 

c  nproj  =  no  of  table  entries 
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c  file  16.2. .. :str(i) ,ptr(i) 

c  str(i)  *  projectile  travel  (cm) 

c  ptr(i)  *  resistive  pressure  (MPa) 

c  interpolate  table  to  find  projectile  resistance  pressure 

c 

c  file  17:  rhli, rkl,rk2 

c  rhli  *  liquid  density  at  zero  pressure  (g/cm**3) 
c  rkl  »  bulk  modulus  at  zero  pressure  (MPa) 
c  rk2  =  derivative  of  bulk  modulus 
c 

c  file  18:  ener,gam 

c  ener  =  chemical  energy  of  propellant  (joules/g) 

c  gam  *  ratio  of  specific  heats 
c 

c  file  19:  sigma ,visk 

c  sigma  *  surface  tension  of  propellant  (dynes/cm) 
c  visk  =  kinematic  viscosity  (cm**2/s) 

c 

c  file  20:  wg,cov 

c  wg  =  molecular  weight  of  the  gases  (g/mole) 

c  cov  =  covolume  of  the  gas  (cm**3/g) 

c 

c  file  21:  pli,pgi 

c  pli  =  initial  liquid  pressure  (MPa) 

c  pgi  *  initial  gas  pressure  (MPa) 
c 

c  file  22:  tdrop  -  droplet  option  switch 
c 

c  if  dropl  then  instantaneous  burning 

c 

c  if  drop2  then  fixed  size  droplets 

c  file  22.1: ddr, pbr 

c  ddr  =  diameter  of  droplets  (cm) 

c  pbr  =  break  pressure  for  change  in  burning  rate  (MPa) 

c  file  22.2:adrl/bdrl 

c  burning  rate  =  adrl  *  (p**dbrl) 

c  if  p<pbr 

c  file  22.3:adr2,bdr2 

c  burning  rate  *  adr2  *  (p**dbr2) 

c  if  p>pbr 

c 

c  if  drop3  then  fixed  size  droplets  wrt  space 

c  droplet  size  changes  with  piston  travel 

c  file  22.1:ddr,pbr 

c  ddr  is  ignored 

c  pbr  *  break  pressure  for  change  in  burning  rate  (MPa) 

c  file  22.2:adrl,bdrl 

c  burning  rate  -  adrl  *  (p**dbrl) 

c  if  p<pbr 

c  file  22.3:adr2,bdr2 

c  burning  rate  =  adr2  *  (p**dbr2) 

c  if  p>pbr 
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c  file  22. 4: ndiam 

c  ndiam  -  no.  of  table  entries 

c  file  22.5... :sfrac(i) , diam(i) 

c  sfrac(i)  =  fraction  of  piston  travel 

c  will  be  normalized  to  max  piston  travel 

c  diam(i)  -  droplet  diameter 

c  the  droplet  diameter  at  any  time  will  be  found  by 

c  interpolating  the  above  table 

c 

c  if  drop4  then  fixed  size  droplets  wrt  space 

c  droplet  size  changes  with  projectile  travel 

c  file  22.1:ddr/pbr 

c  ddr  is  ignored 

c  pbr  ®  break  pressure  for  change  in  burning  rate  (MPa) 

c  file  22.2:adrl,bdrl 

c  burning  rate  *  adrl  *  (p**dbrl) 

c  if  p<pbr 

c  file  22.3:adr2,bdr2 

c  burning  rate  =  adr2  *  (p**dbr2) 

c  if  p>pbr 

c  file  22. 4: ndiam 

c  ndiam  *  no.  of  table  entries 

c  file  22.5. . . :sproj (i) ,diam(i) 

c  sproj(i)  =  projectile  travel 

c  diam(i)  *  droplet  diameter 

c  the  droplet  diameter  at  any  time  will  be  found  by 

c  interpolating  the  above  table 

c 

c  file  23:  tprim  -  primer  option  switch 
c 

c  if  priml  then  primer  instantaneously  goes  to  gas 

c 

c  if  prim2  set  primer  as  liquid  droplets 

c  assume  primer  has  same  properties  as  liquid  propellant 

c  file  23.1:  prprim 

c  prprim  *  pressure  the  primer  should  generate 

c  program  computes  the  primer  mass 

c  set  primer  as  liquid  droplets  in  chamber/ tube 

c  the  initial  gas  pressure  is  assumed  to  have  been 

c  by  primer  that  has  already  burned 

c 

c  if  prim3  then  inject  primer  at  a  constant  rate 

c  assume  primer  has  the  same  properties  as  liquid  propellant, 

c  except  energy  content  may  be  modified, 

c  inject  into  gun  tube, 

c  file  23.1:  rmprim,  tinjec,  fener 

c  rmprim  -  mass  of  primer  (g) 

c  tinjec  *  time  for  complete  injection  (s) 

c  fener  *  factor  multiplied  by  propellant  energy  content, 

c 

c  file  24:theat  -  heat  loss  option 
c 
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c  if  heatl  then  no  heat  loss  to  the  tube  walls 

c 

c  if  heat2  then  convective  heat  loss  to  the  gun  tube 

c  fits  to  viscosity  and  thermal  conductivity, 

c  made  for  HAN1845  from  1500k  to  3500k. 

c  file  24. l:tempw,hlfac 

c  tempw  =  tube  wall  temperature  (K) 

c  hlfac  *  fudge  factor 

c 

c  file  25:  tshock 

c 

c  if  shockl  no  air  shock  term 

c 

c  if  shock2  include  air  shock  effect 

c  file  25.1:airp,airt 

c  airp  =  initial  pressure  in  gun  tube  (MPa) 

c  airt  =  initial  temperature  in  gun  tube  (K) 

c  file  25.2:airgam,airmw 

c  airgam  =  ratio  of  specific  heats  of  gas  in  barrel 

c  airmw  =  molecular  weight  of  gas  in  barrel  (g/mole) 

c 

c  file  26:  ttube 

c 

c  if  tubel  then  standard  lagrange  distribution 

c 

c  file  27:  tinc,htop 

c  tine  *  time  increment  between  output  lines  (s) 
c  htop  *  max  time  step  allowed  in  the  integration  (s) 
c 

c  file  28:  eps,srec 

c  eps  =  error  control  for  time  integration 

c  srec  =  cutoff  between  relative  and  absolute  error  control 

c 

c  file  29:mf,kwrite 

c  mf  =  method  of  integration 

c  usually  mf=22  (stiff  integration  scheme  with  internally 

c  generated  jacobian) 

c  kwrite  =  diagnostic 

c  if  kwrite=l,  then  print  out  after  each  time  step 
c 

c  file  30:tmax 

c  tmax  *  max  run  time  allowed  (s' 

c 

c  file  31:trep 

c 

c  if  repl  then  integrate  once 

c 

c  if  rep2  then  iterate:  change  vent  area 

c  to  adjust  liquid  pressure 

c  file  31.1:  ptar,  vine 

c  ptar  =  desired  liquid  pressure  (MPa) 
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c  vine  =  initial  increment  for  changing  vent  area  (cm* *2) 

c 

c  file  32:tcham 

c 

c  if  chaml  then  compute  chamber  pressure, 

c 

c  if  cham2  then  use  experimental  chamber  pressure, 

c  read  in  from  standard  graphics  file, 

c  file  32.1:  filel 

c  filel  =  file  with  experimental  chamber  pressure 

c  file  32.2:  koll,  kol2 

c  column  koll  =  time  (ms) . 

c  column  kol2  =  chamber  pressure  (MPa) . 

c 

c  if  cham3  then  use  time  derivative  of  chamber  pressure, 

c  integrate  to  get  chamber  pressure, 

c  read  in  from  standard  graphics  file, 

c  file  32.1:  filel 

c  filel  =  file  with  experimental  chamber  pressure 

c  file  32.2:  koll,  kol2 

c  column  koll  =  time  (ms) . 

c  column  kol2  =  derivative  of  chamber  pressure  (MPa/s) . 

c  *********************************************** 
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Below  is  a  listing  of  the  job  stream  for  the  Round  4  model 
(droplet  burning) .  Following  is  a  listing  of  the  output. 

Since  there  are  many  variables  to  keep  track  of,  the  output  is 
on  multiple  files.  Only  the  ones  of  interest  for  a  particular 
problem  need  to  be  printed  out.  Here  only  the  minimal  output  is 
given. 


dr04-30mm  rap-hanl845-l/2  ch-y  hole  block-Drop4-dc=0.95-prim3 


1.88 

252.6 

offset  travel 

3.05 

d4 

114.9 

proj  weight 

2427. 

piston  weight 

84.2 

36.8 

Vl  V2 

43.741 

59.532  0.792 

al  a2  agres 

0.0 

11.498 

avl2  av24 

.2685 

.1080  .1605 

avl  av4a  av4b 

12 

.8687 

nvo  pth 

1.633 

1.93  4.30 

sol  so4a  so4b 

pisl 

piston  resistance 

0.0 

0.0 

1.0 

0.0 

disl 

dis.  coeff.  vs.  pist 

2 

discharge  coeff 

0.0 

0.95 

1.0 

0.95 

fluxl 

fluxl 

projl 

3 


steady  state  bemoulli 
isentropic  flow  into  tube 
projectile  resistance 


0.0 

65.0 

.001 

5.0 

252.6 

5.0 

1.45167 

5600.0  9.45 

rhi  ki  k: 

4287.3 

1.2179 

energy  gam 

71.6 

.06438 

surface  tens: 

23.072 

.707 

mol  wt  gas 

1.03 

.1 

pi  p3 

drop4 

droplets 

0.01 

95.2590 

ddr  pbr 

1.64 

.103 

adrl  bdrl 

1.64 

.103 

adrl  bdrl 

8 

0.0 

.120 

0.1 

.120 

1.0 

.050 

5.0 

.050 

10.0 

.0075 

25.0 

.0075 

45.0 

.0020 

252.6 

.0020 

i  kinematic  viscosity 
covolume 
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0015 


0.70 


prim? 

3.00 

heat2 

300.0 

shock2 


1.0 


0.1  300. 

1.4  28.84 

tubel 

5. 00e-05  1.00e-05 
1.00e-06  1.00e-09 


22  0 
80.0 


repl 

chaml 


inject  primer 
rmprim  tinjec 
heat  loss  to  walls 

tube  temp  fudge  factor 
air  shock 
airp  airt 
airgam  airmw 

standard  lagrange  distribution 
tine  htop 
eps  srec 
mf  kwrite 
tmax 

integrate  once 
compute  p2 
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dr04-30aa  rap-han1845 - 1/2  ch-y  hole  bl ock - D r opt -de *0 . 95 - pr i m3 
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fluxl  steady  state  bernoulli 


rate  of  injection  *  2000. 
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loading  denaity  *  0.90735 
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APPENDIX  C 

JOB  STREAM  FOR  THE  ROUND  25  MODEL 
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Below  is  a  listing  of  the  job  stream  for  the  Round  25  model 
(droplet  burning) .  Following  is  the  summary  page  from  the 
output . 


dr25-30mm 

rap-hanl84 5-1/2 

ch-stag  hole-Drop-dc=0.95-prim3 

1.96 

252.6 

offset  travel 

3.05 

d4 

350.0 

proj  weight 

2427. 

piston  weight 

84.2 

36.8 

vi  v2 

43.741 

59.532  0.792 

al  a2  agres 

0.0 

11.498 

avl2  av24 

.1366 

.0401  .0965 

avl  av4a  av4b 

12 

.8687 

nvo  pth 

1.716 

1.93  4.3.0 

sol  so4a  so4b 

pisl 

O 

piston  resistance 

0.0 

0.0 

1.0 

0.0 

disl 

dis.  coeff.  vs.  piston  travel 

2 

discharge  coeff 

0.0 

0.95 

1.0 

0.95 

fluxl 

steady  state  bernoulli 

fluxl 

isentropic  flow  into  tube 

projl 

projectile  resistance 

j 

0.0 

45.0 

.001 

5.0 

252.6 

5.0 

1.45167 

5600.0  9.45 

rhl  kl  k2 

4287.3 

1.2179 

energy  gamma 

71.6 

.06438 

surface  tension  kinematic  viscosity 

23.072 

.707 

mol  wt  gas  covolume 

1.03 

.  1 

pi  P3 

drop  4 

droplets 

0.01 

95.2590 

ddr  pbr 

1.64 

.103 

adrl  bdrl 

1.64 

7 

.103 

adrl  bdrl 

/ 

0.0 

.090 

0.01 

.090 

0.1 

.025 

1.0 

.020 

4.0 

.0125 

5.0 

.005 

252.6 

.005 

prim3 

inject  primer 

2.50 

.002  0.70 

rmprim  tinjec 

heat2 

heat  loss  to  walls 

300.0 

1.0 

tube  temp  fudge  factor 
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shock2 

0.1  300. 

1.4  28.84 

tubel 

5 . 00e-05  1.00e-05 

1.00e-06  1.00e-09 

22  0 
80.0 

repl 

chaml 


air  shock 
airp  airt 
airgam  airmw 

standard  lagrange  distribution 
tine  htop 
eps  srec 
mf  kwrite 
tmax 

integrate  once 
compute  p2 
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muzzle  vel  (m/s) 
max  v  pis  (m/s) 


913.9 


max  pi  (MPa) 
max  p2  (MPa) 
max  pi  (MPa) 
max  pr  (MPa) 

max  acc  (k-g) 

max  mass  error 
max  energy  error 

ballistic  efficiency 
expansion  ratio  = 

loss  to  tube  walls  = 

run  time  =  6.2 


10.5 

258.6 
192.5 

192.7 

188.7 

37.4 

0.00 

0.29 

=  27.89  % 

14.86 

11.99  % 

nstep  *  1591 
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